X— Systems for Correlation Functions 



J. Caetano a ' 6 ' c , J. Toledo"' 6 



a Perimeter Institute for Theoretical Physics 
Waterloo, Ontario N2L 2Y5, Canada 

Department of Physics and Astronomy & Guelph- Waterloo Physics Institute, 
University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 

c Centro de Fisica do Porto e Departamento de Fisica e Astronomia, 
Faculdade de Ciencias da Universidade do Porto, 
Rua do Campo Alegre, 687, 4169-007 Porto, Portugal 

jd. caetano . s , Jonathan. campbell . toledo@gmail . com 



Abstract 

We consider the strong coupling limit of 4-point functions of heavy operators in M = 4 SYM 
dual to strings with no spin in AdS. We restrict our discussion for operators inserted on a line. 
The string computation factorizes into a state-dependent sphere part and a universal AdS 
contribution which depends only on the dimensions of the operators and the cross ratios. We 
use the integrability of the AdS string equations to compute the AdS part for operators of 
arbitrary conformal dimensions. The solution takes the form of TBA-like integral equations 
with the minimal AdS string-action computed by a corresponding free-energy-like functional. 
These TBA-like equations stem from a peculiar system of functional equations which we call 
a x-syst em - 111 principle one could use the same method to solve for the AdS contribution 
in the Appoint function. An interesting feature of the solution is that it encodes multiple 
string configurations corresponding to different classical saddle-points. The discrete data 
that parameterizes these solutions enters through the analog of the chemical-potentials in 
the TBA-like equations. Finally, for operators dual to strings spinning in the same equator 
in S 5 (i.e. BPS operators of the same type) the sphere part is simple to compute. In this 
case (which is generically neither extremal nor protected) we can construct the complete, 
strong-coupling 4-point function. 
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1 Introduction 

One of the most interesting objects to study in an interacting quantum field theory is the 
four-point function. In conformal field theories these are the first iV-point functions whose 
spacetime dependence is not explicitly fixed by conformal symmetry. The computation of 
this correlator is generally highly nontrivial and obtaining explicit expressions outside the 



3 



perturbative regime is typically impossible. However, the advent of the AdS/CFT corre- 
spondence has made it possible to access the strong coupling limits of special QFT's. In 
particular the correspondence [1] maps the strong coupling limit of M = 4 SYM to a theory 
of classical strings moving in an AdS^ x S 5 background. This allows the computation of 
leading strong coupling behavior of certain 4-point functions in M = 4 SYM, which is the 
main purpose of this paper. While we have mostly focused on the 4-point computation, we 
note that the method used below could in principle be applied for the corresponding iV-point 
computation. 

There are also interesting indirect applications. One of our main motivations for this 
work is the possibility of exploring the spectrum and the structure constants of Af = 4 SYM 
through the operator product expansion (OPE) of the four point function. Another poten- 
tially interesting application of this calculation is in the weak-strong coupling connections 
reported in [2, 3, 4, 5, 6] that suggests the existence of an underlying common structure in 
both regimes. 

At strong coupling, the problem of computing the correlation function is that of finding 
the area of the minimal surface in AdS§ x S 5 that goes to the AdS boundary at the operator 
insertion points x a . In this paper, we compute the AdS part of the correlation function for 
arbitrary heavy scalar operators inserted along a line. The method used here is inspired by 
the integrability techniques originally developed for the Null Polygonal Wilson loop problem 
[10] and later applied to the computation of three-point functions [8, 9, 11]. As in these 
previous applications, integrability allows one to compute the minimal AdS action without 
knowing the explicit classical solution. For the four point correlation function the connection 
with Hitchin systems and the formalism developed in [12] is used intensely. As in [10, 8, 9] 
the starting point of the method is the map of the string equations of motion in AdS to a 
certain the auxilliary linear problem by Pohlmeyer reduction. Ultimately the solution takes 
the form of a set of functional equations that we call a x-system. These functional equations 
are similar in spirit to the F-system appearing in [10] and which naturally arise in the solu- 
tions of integrable QFT's. 

For some specific BPS operators dual to strings spinning on the same great circle of S 5 
the sphere contribution is well know. In this case we can construct the full strong-coupling 
correlation function. We emphasize that these 4-point functions are generically neither ex- 
tremal nor protected. Complete, non-protected results for correlation functions of heavy 
operators at strong coupling are quite rare. For example, in [8], the AdS part of the three- 
point function is computed, but the only case for which the sphere is known (BPS operators) 
is protected. The only complete, non-protected result that we know of is [11] where the 
strong-coupling three-point functions of GKP strings is computed. 1 

The layout of the paper is as follows. In the section 2 we start by giving the general 
strategy of this work and discuss some physically relevant aspects of the problem. Then, in 
section 3 we write the AdS part of the correlation function in terms of objects which are 
naturally computed using the integrability of the string equations of motion. In section 4 we 
introduce a formalism that will lead to the x-system which is the set of functional equations 

1 Using the results of this paper it may be possible to extend the results of [11] to the complete iV-point 
functions of GKP strings at strong coupling since the mathematical problem is similar to the one treated 
here. 
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that allows one to compute the minimal AdS action. In section 5 we compute explicitly the 
AdS part of the correlation function using the x-system and explain the mechanism by which 
the dependence on the spacetime points emerges. We also present some numerical tests of 
the results. In section 6.1, we compute the sphere part for the specific case of BPS operators 
with large charges. In section 6.2, we discuss the saddle points of the fourth insertion point 
of the correlation function and perform some numerical studies on this issue. In the same 
section, we study the extremal limit of the correlation function which is an analytical test of 
the result. In section 7 we conclude and discuss some open problems. 

Finally, because the computation involves many intermediate steps in the flow-chart below 
we summarize for the reader the basic components of the method. The flow-chart provides 
a map of the paper and summarizes the main steps of the computation which are shown in 
the individual boxes. Each box contains a reference to the relevant section where that part 
of the computation is discussed. 
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2 Four point function generalities 



For large 't Hooft coupling A, the semi-classical computation of correlation functions corre- 
sponds to the evaluation of the AdS§ and S 5 actions for classical solutions with the topology 
of a four punctured sphere. The boundary conditions are that the solution close to each 
puncture P a , which is associated with the gauge theory operator O a (x a ), approaches the 
AdS boundary at the point x a in the same way as a 2-point function involving O a (x a ) and 
some other heavy, scalar operator. In this paper, we study the simplest case where the op- 
erators are inserted on a line in IR 4 . This implies that the string solution is contained in a 
Euclidean AdS2 subspace of AdS$. Moreover, there is only one independent cross-ratio. 
The conformal symmetry of TV = 4 constrains the four-point correlation function to take the 
form 

4 

(0 1 (x 1 )0 2 (x 2 )0 3 (x 3 )0,(x 4 )) = f(u) n(^) Aab > (1) 

a>b 

where x a b = x a — Xf,, A a is the dimension of operator O a , A ab = (J2 C /3 — A a — A b and 
u is the conformal cross-ratio 

^14^23 fct \ 

u = (2) 

£12X34 

Both the AdS and sphere contributions contain divergences as the string approaches the 
position of the operators at the boundary of AdS, which requires a cut-off z — S. To describe 
the world-sheet we use complex variables w,w. On the 4-punctured sphere, the physical cut- 
off £ corresponds to cutting out small disks of radius e a around each puncture P a at w a . 
Ultimately, we will need to establish a precise relation between the cut-off's e a and £. As we 
will review later, this is possible given the data accessible from integrability [8]. 

In this paper, we will consider operators with charges scaling as a/A, and without spin 
in AdS. Following the prescription developed in [7, 8], we account for the states in the 
sphere by introducing an extra contribution of wave-functions. Therefore, the semi-classical 
four-point function is given schematically by 

/^^W^e-^W^™ (3) 

where the actions are evaluated on a classical (Euclidean) string solution approaching the 
boundary of AdS at the positions of the insertion points x a . 

In principle, there is an integral over all four insertion-points on the worldsheet. In (3) 
we only integrate over the insertion W4 since the position of the other punctures can be fixed 
by conformal transformations. Since we are considering the A — > 00 limit, one can evaluate 
the integral over W4 by saddle point and the end result is the integrand of (3) evaluated at 
the dominant saddle point. 

Let us consider the issue of the saddle point in some detail since it will be an interesting 
aspect of our computation. There are two issues here: the positions of the operators on the 
boundary and the positions of the insertion points on the sphere. We can use the target- 
space conformal symmetry to place three of the operators at x% = 1, x 2 = 00, x 3 = — 1 and 
the world-sheet conformal symmetry to fix Wi = 1, w 2 = 00, w 3 = — 1. The position x 4 is an 
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X4 X3 X\ X2 X4 X3 X\ X2 




Figure 1: Insertions on the 4-punctured sphere. The gray ball represents the world-sheet (the 
complex plane plus the point at infinity, or simply 'the sphere') and the black boundary of the ball 
represents the equator of the sphere. The points w a are the punctures on the sphere corresponding 
to the operators inserted at the positions x a at the boundary of AdS2, which is represented by the 
strait line. We fix the points wi, W2, 103 and x±, x%, X3 using the world-sheet and target-space 
conformal symmetry respectively. The position of the fourth insertion W4 should be fixed at the 
dominant saddle point w\ of the integrand of (3). By symmetry we expect this saddle point to also 
be along the real axis, and thus we have a notion of an ordering of the 4 punctures (but see footnote 
2). In particular, there is three distinct ranges for the location of w\. Consider the ordering of the 
x a shown in this figure. If the dominant saddle point is located between W2 and W3 (as in panel 
A) then the insertions will not cross and the string embedding will look schematically like the one 
shown in figure 2A. If the dominant saddle-point is located between W3 and wi (as in panel B) 
or between w\ and u>2 then the insertions cross each-other and we expect the string embedding to 
look like the one shown in figure 2B. 



input since we can put O4 anywhere along the line that contains ^1,2,3- On the other hand, 
once we choose X4 the position of the fourth puncture is fixed at W4 = w* A by the saddle-point 
condition. By symmetry we expect the dominant saddle-point to be located on the real axis 
and in this case we have a notion of an ordering of the punctures. 2 In particular, there are 
three possible in-equivalent orderings depending on the position of W4. Figure 1 shows two 
of these possibilities. If the ordering of the x a is the same as the w a then the insertions do 
not cross each other, as in figure 1A. If the ordering of the x a is different from that of the 
w a , then the insertions will cross as in figure IB. These two possibilities lead to two types 
of string embeddings with distinctly different properties as is shown in figure 2. We will see 
that two types of solutions arise naturally in our construction. We are able to characterize 
the qualitative features of the spacetime embeddings and compute the minimal AdS action 
of both types of solutions. We will return to this topic below. 

2 We have confirmed numerically that there are saddle points along the real axis. There may also be 
the possibility saddle points off of the real axis and occurring in complex-conjugate pairs, but we have not 
investigated this possibility. This issue certainly requires further investigation. 
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X4 X3 X\ X2 £4 X3 X\ X2 X4 X3 X\ X2 




A B 

Figure 2: Two different possible string embeddings in AdS 2 which obey the required boundary 
conditions. These two solutions are shown in panels A and B. The center panel shows how to 
generate the configuration of panel B from that of panel A by interchanging the order in which 
the insertions on the sphere attach to the boundary; this interchange results in the characteristic 
folding shown in the embedding of panel B. These two types of solutions arise from the possibility 
that for a given choice of operator insertion points x a the insertion point w\ (see figure 1) can be 
located in any of the three possible intervals (w2,ws), (^3, w\), (wi,W2)- 

3 AdS2 Pohlmeyer reduction 

In this section we briefly review the Pohlmeyer-reduction process. We begin with a discussion 
of the string equations of motion and the stress-energy tensor, which is the starting point 
of the reduction. We then introduce the function 7 in terms of which the AdS Lagrangian 
can be written. It turns out that 7 satisfies a non-linear but scalar equation of motion 
that is a modified version of the well-know sinh-Gordon equation. Next we show how the 
different types of string embeddings discussed in section 2 are encoded though the boundary 
conditions imposed on 7. Finally we use the function 7 to write the AdS action in a form 
where integrability is more readily applied. 

3.1 Equations of motion and stress-energy tensor 

Recall that we can consider (euclidean) AdS2 as a surface in R 1,2 obeying the constraint 

Y-Y = (Y 1 ) 2 -(Y 2 f + (Y 3 ) 2 = -l. (4) 
We write the action for a string in AdS 2 as 

S=^Jd 2 a [d a Y -d a Y + \{Y -Y + 1)] (5) 

and the resulting equations of motion as 

UY = (BY ■ BY) Y (6) 

The first term in the action is just the free string action in R 1 ' 2 ; the second term is a Lagrange 
multiplier term that imposes (4). 
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X 

X 

Figure 3: Schematic analytic structure of T. The blue dots represent the (double) poles of T 
at locations w a and corresponding to the operator insertions O a (x a ). The yellow crosses indicated 
zeros of T. We have fixed the positions of wi, u>2 and W3 using the world-sheet conformal symmetry. 
We have arbitrarily placed 1114 in the interval (ws,wi) although generically the saddle-point w\ can 
be located in any of the three intervals along the real axis. 



The equations of motion (6) must be supplemented by the Virasoro constraints and 
boundary conditions. The Virasoro constraint requires T^ds + T s = 0. In particular, the 
AdS contribution to the stress-energy tensor does not vanish. Fortunately the boundary 
conditions allow us to completely fix the form of Tacls — —Ts- Here we are interested 
in solutions with the topology of a four-punctured sphere where the punctures are at the 
position of the operator insertions and thus the boundary conditions give the behavior of 
the string solutions near the insertion points. The correct prescription is to demand that 
the string goes to the boundary at the insertion points. Furthermore, it should approach 
the boundary in a specific way as dictated by the vertex operators. The behavior of the 
solution near the boundary will be dominated by the operator inserted there, independent 
of the properties or number of other operators inserted at different points. This means that 
the behavior near the insertion points can be determined from the 2-point function, where 
the string solution is know explicitly. From the explicit solution for the 2-point function one 
finds that the desired property of the solution near insertion point w a is [8] 

{dYf = I» ~ Al (w w a ) (7) 

where T (w) is the holomorphic component of T^ds- The corresponding property also is 
required for the anti-holomorphic component T(w). Thus we know that T should be an 
analytic function on the (4-punctured) Riemann sphere with double-pole singularities at the 
punctures. This fixes T to be a specific rational function. 

First consider the denominator of the rational function T. The polynomial in the denom- 
inator is determined by the positions of the insertions. Three of the insertions can be fixed 
by conformal symmetry, leaving one final insertion. The integrand of (3) will be a function 
of this final insertion point. In the limit \/A — > 00 the integral localizes at the saddle point 
w± = w\, thus fixing completely the denominator of T. 

Now consider the numerator of T. Without loss of generality we can consider the case 
where there is no insertion at infinity since we can perform a transformation that maps any 
arbitrary point to infinity. Then the polynomial in the numerator can be at most of degree 4 
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(otherwise T would not be regular at infinity) and therefore it is characterized by 5 param- 
eters. Four of these parameters are fixed by the condition (7). The final unfixed parameter, 
which we will call U, parameterizes the single cross-ratio of the four operators (recall that 
four points in a line have only one independent cross-ratio). The precise map between the 
parameter U and the cross-ratio u is quite involved but fortunately we will not need it since 
the cross-ratio will be encoded in the x-system in a simple way. The analytic structure of T 
is shown schematically in figure 3. We will use this sort of figure to represent T throughout 
this paper. 



3.2 The function 7 

Our objective is to evaluate the AdS part of the string action. In Poincare coordinates the 
on-shell action becomes 3 

„, , „, , dxdx + dzdz r—= , . , 

dY ■ dY = = V TT cosh 7 (9) 

z 2 

where the above formula defines the function j(w, w). It follows from the equations of motion 
that 7 satisfies the modified sinh-Gordon equation 

dd-f = VTfsinh'j. (10) 

It is well known that this equation is classically integrable, and in what follows we exploit 
this integrability to compute the AdS action. 

Now let us determine what boundary conditions should be imposed on 7. For the 2-point 
function 7 = 0. Recall that the string solution should approach that of the 2-point function 
as the string approaches the boundary at the operator insertion points x a . Therefore we 
should require that 7 — > as w — > w a [8]. Furthermore, in order to have a non-singular 
world-sheet metric the right-hand side of (9) should never vanish. Thus when T has a zero 
7 must have a logarithmic singularity to cancel it. In summary, the boundary conditions on 
7 are 

7 ->• ±^logTf (w^z a ) (11) 

7^0 {w->w a ) (12) 

where z a denotes a zero of T and w a a pole of T. Notice that the regularity of the world- 
sheet metric does not fix the sign of the logarithmic 'spike' in (11) and, in principle, different 
choices are possible at each zero (recall that generically T will have 4 zeros for the 4-point 
function, as follows from the discussion of the previous section). These different choices 
correspond to different string solutions having differing properties, and generically different 



3 The AdS 2 Poincare coordinates are given by 



Y i = -±(l-x 2 -z 2 ), = 1 (l + ^ + z 2 ), F 3 = -. (8) 

2z x ' 2z x ' z 
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Figure 4: Contours where 7 = based on the choice of signs in equation (11). These contours are 
shown schematically by the black curves. The label u (d) at a zero indicates the choice of sign + 
(— ) in equation (11). We give a detailed discussion of why these are the only possible structures 
for these contours in appendix C. The key in relating these figures to the embeddings in figure 2 is 
that contours on the world-sheet where 7 = map onto folds of the embedding. 



total action. We will refer to the spikes with the + (— ) sign as w-spikes (d-spikes). We will 
now describe how the choice of these signs is related to the string embeddings shown in figure 
2. 



3.3 Spikes, fold-lines and string embeddings 

As mentioned in the previous section there are 4 zeros of T and at each zero we have a Z2 
ambiguity (see equation (11)) in the choice of spikes of 7. A priori there are 2 4 different 
choices for the spikes. However, it turns out that there are only 2 distinct choices that 
correspond to target-space solutions with the desired properties. These two possibilities are 
shown in figure 4. A discussion of why these are the only two possible choices is given 
in appendix C. 4 These two different possibilities correspond precisely to the two different 
possible string solutions shown in figure 2. The key ingredient in making this correspondence 
is the observation that contours on the world-sheet where 7 = correspond to fold-lines in 
the string embedding (see appendix C). The location of these contours is directly connected 
with the choice of spikes. For example, between a w-spike and a <i-spike we know that there 
must be at least one such contour. In figure 4 the 7 = contours are indicated by the black 
curves. In appendix C we discuss in detail how the structure of these contours is inferred 

4 Thc main ideas are: first, configurations related by 7 — > —7 are not distinct since this is a symmetry of 
(10), and second, one should choose the spikes such that 7 — > —7 under reflection about the real axis. See 
appendix C. 
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from the choice of spikes. 

Let us describe in more detail how we relate the two spike configurations in figure 4 to the 
two string embeddings in figure 2. As mentioned above, the key ingredient is to study the 
fold lines in the two figures. First consider the target-space solution. In figure 2A it is clear 
that there is a single fold-line that runs through each of the punctures in sequence. That is, 
there is a fold-line directly connecting 24 with 23 then 23 with xi, etc. This is in agreement 
with the fold-structure implied by figure 4A since for this choice of spikes we can deduce 
that there is a single contour where 7 = running along the real axis connecting to W3 
then u> 3 to wi, etc. Thus the spike configuration of figure 4 describes a string embedding of 
the type shown in figure 2A. 

Now consider the folds of the embedding shown in figure 2B. Insertions X\ and 23 are 
both connected by fold lines directly to the insertions 22 and 24. Furthermore, 22 and 24 are 
connected to each other by two fold-lines. This is because this configuration is double-folded 
along that line, as one can see from the construction shown in the center panel of figure 2. 
All of this is in perfect agreement with the fold-structure implied by figure 4B. In particular, 
for this choice of spikes both W\ and w 3 are directly connected to w 2 and w% by contours 
where 7 = 0. Moreover, w 2 and w* A are connected by two contours where 7 = 0, precisely 
corresponding to the double-fold line connecting 22 and 24 in figure 2B. 

Let us comment on a subtle point regarding figure 4. Note that we have placed the saddle 
point w\ in different intervals in the two figures. On one hand, we should do this in order to 
be in agreement with figures 1 and 2. However, as we will see, given a cross-ratio u and the 
saddle point wl, the fold structure is fixed. So, to compare two different fold structures for 
a given cross ratio we are forced to place the saddle point wl in different intervals. This is 
in perfect agreement with the intuitive perspective of figures 1 and 2. We will return to this 
point in section 6.3. 



3.4 The action as a wedge product 

We will now return to the computation of the minimal AdS action (see equation (3)). Ex- 
plicitly, the quantity we want to evaluate is 

f dxdx + dzdz , . 

X JZ\{e a } z 

where £\{e a } denotes the sphere with small disks of radius e a cut out at each puncture. These 
cut-offs are not independent and are all fixed in terms of the single target-space cut-off z = £; 
this is important in recovering the spacetime dependence of the correlation function and we 
will return to this point below [8]. It is convenient to separate the action into a piece that is 
independent of the cut-offs, and a piece where the dependence can be explicitly evaluated. 
This can be done because the solution near the punctures is know. In particular, we may 
write [10, 8] 



A= -— [ v / TT(cosh 7 - 1) - — / Vrf (14) 

To extend the integration to the full sphere in the first term we have used the fact that 
the action (9) goes like y/TT near the punctures as follows from (12). We will refer to the 
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first and second term in (14) as Af in and A div respectively. Since T is known A div can be 
evaluated explicitly in terms of the e , but to eliminate the e a in terms of £ requires detailed 
information about the string solution itself. Fortunately, the tools necessary for computing 
Afin will also provide the necessary information to complete the calculation of A^ v . Thus, 
let us focus for the time being on the computation of Af in and return to A^ afterwards. 

We would like to write Af in in a form where the integrability of (10) is more readily 
usable. Following [10, 8] we introduce the forms 

to = Vfdw (15) 

1 r= 11 

7] = - v T (cosh 7 — 1) dw + -— = (d^) dw (16) 

2 4 y/T 

and then from a direct computation it follows that 



A fin = - luAr] (17) 
J s 

where £ denotes the double cover of the sphere defined by y 2 = T(w). Extending the 
integration from E to E simply involves a factor of 2 since each form is odd under sheet- 
exchange. An important property of these forms is that they are both closed. The form u 
is clearly closed since it is holomorphic, and the closure of rj follows from the equations of 
motion for 7. Notice that (17) would be true for any choice of the dw component of 77. The 
specific coefficient appearing in (16) is necessary for the closure of the form. 

Now we would like to apply the Riemann bilinear identity (RBI) to reduce (17) to one- 
dimensional integrals over cycles on S. There are two caveats in doing this - the singularities 
in u and the singularities in 77. These issues where resolved in [8], and we follow the approach 
used there (see [8] for a more detailed treatment and also [9] for a different approach). The 
basic idea of the RBI is to write one of the forms of the wedge product as an exact form, 
uj = dF where F = J Po 00 , which is always possible on a Riemann surface minus some contour, 
L. In the present case u has single poles and thus F will have logarithmic cuts which need to 
be accounted for. A way to side-step this complication is to spread the single poles in u into 
a small square-root cuts such that F has only square-root cuts and no singularities. The cost 
of doing this is that the genus of £ increases, but the upside is that the application of the 
RBI is simplified. This takes care of the singularites in u. Now consider the form 77 which 
behaves as 77 ~ (w — z a )~ 5 / 2 near the zeros of T. The prescription of [8] is to remove the 
points z a from the domain by taking L to be the sum of the standard contour for a Riemann 
surface of genus g and small contours C a encircling the points z a . The integrand of (17) can 
then be written as d (Frj) (since drj = on the domain) and then Stokes theorem can be 
used to reduce the surface integral to a line integral over the usual boundary of the genus g 
Riemann surface and the contours C a . The end result is that each boundary C a contributes 
a correction of 7r/12 to Afi n in (17) while the integral over the boundary of £ gives the usual 
sum over cycles on S and thus we have [8] 

A fin = (number of zeros) ^ - l -(^ju^j l£ \ J> r^j (18) 
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where {7 a } is a complete basis of cycles on E and I a b is their intersection matrix. For the 
four-point function there is generically 4 zeros and 4 poles. When we spread the four poles we 
introduce an additional 4 cuts and thus the surface is genus 5 and there will be 5 a-cycles and 
5 b-cycles; that is {^ a } = {7 ai , 7^ , 7 a2 , j as , 7& s }. The main point is that we have reduced 
the computation of the surface integral (17) into a sum of 1-dimensional cycle integrals of 
a closed form. Such integrals are precisely what integrability is good at computing. In the 
following section we will see how to compute the cycles § i] by exploiting the integrability 
of (10). 

4 The linear problem 

To compute the ^-cycles appearing in (18) it is useful to consider the linear problem associ- 
ated with equation (10). Consider a function ip obeying 

(9+J tt )^ = 0, («9 + J„)V = (19) 

where the components of the connection J = J w dw + Ju,dw are given by 

J w = A w + J iB = A €! + (20) 

where A and $ are independent of the spectral parameter £ and given in terms of 7 and T, 
T. We give the explicit forms of A and $ in appendix A. Note that we will frequently write 
the spectral parameter as £ = e 9 . 

Compatibility of equations (19) for all £ is equivalent to the flatness of J, which is satisfied 
if 7 obeys the equation of motion (10) and T (T) is purely holomorphic (anti-holomorphic). 
In the following section we will discuss the relation between the solutions of the (19) and the 
77-cycles appearing in (18). 

4.1 Basic properties 

There are a few aspects of the linear problem which will be essential for the following analysis. 
Let us comment on each of them in turn. 

• Solutions near punctures. Using (7) and (12) one can show that near the punctures P a 
there are two linear-independent solutions of the form (see Appendix A) 

= (T/f) 1/8 e^f w ^ + ^\±) (21) 

~ (W-Waf^ 1 -* (w - Waf^+l |±) ( 2 2) 

where |±) are the eigenvectors of a 3 . Notice that there is a solution that is exponentially 
big and one that is exponentially small as one approaches the puncture P a . 5 

5 In going from (21) and (22) we have been careless about the branches of ui. In particular, we may choose 
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• 'Small' solutions. Demanding that a function is both a solution of the linear prob- 
lem and also small at some puncture P uniquely defines that solution (up to overall 
normalization). Thus there is a family of 'small' solutions s a each of which is small 
at puncture P a . On the other hand, specifying that a solution has the big asymptotic 
near P does not uniquely determine the solution since one could create another solution 
obeying the same boundary conditions by adding an arbitrary multiple of sp. 

• 7j2 symmetry and 'big' solutions. Even though one cannot uniquely specify a solution 
by demanding that it has the big asymptotic near P, there is nevertheless a special 
solution big near P that is uniquely defined. This follows from the Z 2 symmetry of the 
connection (20) which is given by 

UJiOU- 1 = J(e^) (23) 

where U = ia 3 . This symmetry implies that if sp(£) (we are suppressing the w, w 
dependence) is the solution to (105) small at P then 

s P = a 3 s P (e~ in (24) 

is another solution of the linear problem. Moreover, from (21) it follows that §p is big 
at P. Thus we have a second uniquely defined family of solutions s a , each of which is 
big at puncture P a . 

• Products of solutions. Given two solutions of the linear problem ipx and ip2, there is a 
natural SL2 invariant inner product 

(ViAV 2 ) = Det[{V>i,V2}] (25) 

This inner product is equivalent to the Wronskian of the two solutions. Important 
properties of this Wronskian are that it is independent of w and w, and thus only 
depends on the spectral parameter £. Further, the product will vanish if the two 
solutions are linearly dependent. 

Now that we have introduced these basic facts of the linear problem, we can state what is 
perhaps the key ingredient in the whole computation. 6 We claim that the £ — > expansion 
of the inner product of two small solutions is the following[10, 8] 

(s a A s b ) ~ exp 



a particular branch at some P a such that the near-puncture solutions take the form (22) but then if we 
smoothly continue yT to some other puncture Pf, it is possible that the small and big solutions correspond 
to the opposite components from the small and big solutions at P a . This will be very important below, since 
it will usually be the case in the construction we will use. 

6 To our knowledge the following fact first appear in [10]. Later it was used in [8, 9] for 3-point function 
computations. We give a derivation in appendix B; we thank Pedro Vieira and Amit Sever for explaining 
the basic components of the derivation used for [10]. A different derivation appears in [9]. 



(26) 
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where r\ is precisely the same form (16) that appears in the action formula (18) and zu ab , 
w a b are explicitly known in terms of integrals of u and u. 7 A derivation of (26) is given 
in appendix B. The point is that by computing the inner products (s a Asb)(£) we can 
extract the "puncture-puncture" integrals r] by extracting the O (£) coefficient of this inner 
product. All of the ^-cycles appearing in (18) can be written in terms of linear combinations 
of these puncture-puncture integrals. Thus, we can compute area (18) by computing the 
inner products (s a A s&). The rest of this section is devoted to explaining how we compute 
such inner products using techniques from integrability. 

4.2 Defining solutions globally 

Let us now comment on how to globally define the small solutions. Suppose that we want 
to construct the small solution sp away from puncture P, say at some generic point A. 
To do this we need to use the connection to transport the solution along some path from 
the neighborhood of P to the point A. However, it is clear from (21) that the solutions 
of the linear problem have non-trivial monodromies around the punctures and therefore 
homotopically different paths on the 4-punctured sphere will result in different values of the 
small solution at A. In other words, solutions of the linear problem live on a (generically 
infinite-sheeted) Riemann surface with branch points at the punctures. For the purposes of 
calculating it is convenient to fix some conventions for dealing with the multivaluedness of 
the solutions. We first define the sheets by cutting the Riemann surface as shown in figure 
5. The cuts all join at a common point and the monodromy about that point is the identity 
since a path passing through all the cuts is contractable on the sphere. We then define the 
value of the small solution associated with puncture P at some point A as follows. Draw 
any curve from the neighborhood of P to A. In the neighborhood of P one starts with Sp. 
For every time the path crosses a cut emanating from some puncture Q in the clockwise 
(counterclockwise) sense attach a factor Mq (Mq 1 ). 8 In this way, if we transport along a 
path that is homotopically equivalent to a path that does not cross any cuts then the small 
solution at A will be sp\a- If the path crosses the cut emanating from puncture Q once in 
the clockwise sense, then the value of the small solution at A will be (MqSp)\a, and so on 
(see figure 5). In the case when sp is transported around the puncture P one can see from 
the explicit form (21) of sp near P that the result will be multiplication by a constant. That 
is 

MpSp = fipSp (27) 
MpSp = jlpsp (28) 

7 To be more precise, this expansion will be true for certain s a and s& depending on certain conditions 
stemming from the form of T and also depending on the value of Arg(£). Furthermore, the contour of 
integration will be precisely defined by these conditions. We will discuss these conditions in detail below. 

8 Note that the result of a monodromy can be expressed as the linear map M since both s and Ms are 
solutions of the linear problem. Therefore they can both be expanded in terms of two linearly independent 
solutions of the linear problem, and thus they are related to each other simply by a linear map, or in other 
words simply by multiplication by some matrix, M. 
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(M 3 s 4 )U = (M 2 - 1 M 1 - 1 M 4 - 1 s 4 )U 




Figure 5: Our conventions for denning the solutions globally. The dashed blue lines emanating 
from the punctures indicate the conventions for 'cutting' the full Riemann surface, thus defining 
the sheets. The red lines indicate the parallel transport of a solution from P4 to the point A 
along three paths. Two of the paths are homotopically equivalent due to the triviality of the total 
monodromy M4M3M2M1 = 1 (which follows from the fact that any path encircling all the punctures 
with the same orientation is contractable on the sphere). The third path is homotopically distinct 
from the other two, and thus the value of the solution at A will differ by monodromy factors. 



so that sp and 5p are eigenvectors of M P with eigenvalues \ip and jlp = l/fip respectively. 
One cannot repeat such an analysis to evaluate MqSp since genericaliy one does not know 
the explicit form of sp in the neighborhood of Q. 

4.3 WKB approximation and WKB Curves 

As we will discuss shortly, it will be essential to have control over the £ — > 0, 00 asymptotics of 
the inner products (sp A sq) (£). It is clear from (19) — (20) that these are both singular limits, 
and the basic idea of extracting this singularity - which is called the WKB approximation - is 
as follows. 9 As discussed above, we have good control over the solutions in the neighborhood 
of the punctures. Thus we would like to study, in the limits £ — > 0, 00, the transport of small 
sp along a curve w(t) which connects a neighborhood of a puncture P with a neighborhood 
of another puncture Q. Let us consider the transport away from P (see figure 6). We will 
discuss the £ — > limit since the £ — > 00 limit is similar. 

At any point in £ the matrix $ has the two eigenvalues =Fu;/2 = dw (which are 

single valued on the double cover £), and thus we can choose a gauge along w(t) where $ is 
diagonal and given by 

_ I f-uj 0\ _ 1 (-Vfdw 

9 See appendix B or [12] for a more detailed treatment. 
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Figure 6: Transporting sp away from P along w(t). We have chosen the branch of $ in (29) 
such that sp oc |+) near wp. In other words, we have chosen the branch of <1> such that 

Re((+| — <& w / £dw\+)) = Re (dw\/T/^\ > for duu pointing along w(t) away from P and thus 
exp ( f™y^ dwy/T/^\ is exponentially diverging as £ — > 0. 



In the limit £ — > some component of $/£ will dominate and thus the leading contribution 
to s a at some point w along w(t) will be given by 

e-M^la) (30) 

where the value of a — ± depends on the branch of $ we have chosen (recall that |±) are 
the eigenvectors of a 3 ). This is the singular contribution in the limit £ — > for the same 
reason that it is the small solution - namely, because 

Re (H (-$/£) |a»>0 (31) 

along a path traveling away from P a . The basic statement of the WKB approximation is 
that so long as we transport along paths such that (31) is true along the whole path then 
the leading contribution to sp in the £ — )■ limit is indeed given by (30). In other words, as 
long as we transport along curves satisfying (31) everywhere, then we can reliably extract 
the singularity as £ — > as it is simply given by (30). The curves along which (31) is satisfied 
most strongly are those for which 

Im(H (-$/£) k)) = (32) 

Curves satisfying this condition are called WKB curves. If we transport along some curve 
satisfying (32) for Arg(£) = 0, then the condition (31) will be satisfied for Arg(£) e 
(0 — 7r/2, + 7r/2). In fact, we will need to control the asymptotics of sp in precisely such 
a wedge of the £-plane, and thus we should always transport along WKB lines. We will 
give the a very brief overview of the properties of these lines in the next subsection. For a 
detailed treatment see [12]. 
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A: generic point 



B: pole 



C: zero 



Figure 7: Local structure of WKB lines in the neighborhood of, A: a generic point; B: a double 
pole of T; C: a simple zero of T. In the case of a generic point the WKB curves form continuous 
non- intersecting lines. In the case of a singular point they form logarithmic spirals for generic values 
of Arg(£). The exact nature of these spirals will not be important. What is important is that the 
singular points act as sources/sinks of WKB curves. In the case of a zero, there are three special 
WKB curves that asymptote to the zero which are the red curves in panel C. These special curves, 
called separating curves, determine the global structure of the WKB foliation. 



4.4 WKB triangulation 

As we discussed in section 4.2 we define the solutions of (19) globally by transporting along 
specific paths. Transport of solutions along homotopically equivalent paths will lead to the 
same result, whereas transport along homotopically inequivalent paths generically will give 
different results. For this reason it is useful to set up a system of fiducial paths between the 
punctures which we will use to globally define the solutions. Because we will need to control 
the large/small £ asymptotics of the Wronskians, it is best to choose these paths to be WKB 
curves - i.e. curves satisfying (32). 

We will first consider the local structure of WKB curves. In the neighborhood of a generic 
point on the punctured sphere the WKB curves are smooth and non-intersecting (see figure 
7A). In the neighborhood of a (double) pole of T the WKB curves follow logarithmic spirals 
that asymptote to the singular point (see figure 7B). All that will be important here is that 
the poles act as sources / sinks of WKB curves but the exact nature of these spirals will not 
be important. Finally, working in the neighborhood of a simple zero of T one can show that 
there are three special WKB curves that asymptote to the zero and which govern the WKB 
lines near the zero (see figure 7C). 

Now consider the global structure of the WKB curves. All WKB curves fall into one of 
the following types [12] 

• Generic WKB curves which are those that asymptote in both directions to a pole of 
T (potentially the same one); 

• Separating WKB curves which asymptote to a zero of T in one direction and to a pole 
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Figure 8: Global WKB structure for an example with 4 punctures. The separating curves are shown 
in black. In one cell we show several examples of homotopically equivalent curves (shown in gray) 
that sweep the cell. Each cell defined by the separating curves has a 1-parameter family of such 
curves. By choosing a representative curve from each family we obtain the triangulation shown in 
figure 9. Notice that near each puncture (the blue dots) we see the spiral structure shown in panel 
B of figure 7 and near each zero (yellow x ) we see the local structure shown in panel C of figure 7. 



of T in the other; 

• Finite WKB curves which are closed or asymptote in both directions to a zero of T 
(potentially the same one). 

We will now describe how we use the WKB curves to set up a system of fiducial paths, or 
more specifically, a triangulation. By triangulation we mean a triangulation of the punctured 
sphere with all vertices at the punctures and at least one edge incident on each vertex. 
Consider fixed T and Arg(£) such that there are no finite WKB curves (this can always be 
done since such curves only appear at special, discrete values of Arg(£)). First draw all of the 
separating WKB curves - there will be 3Nz of these, where Nz is the number of zeros of T 
(since for the moment we are not allowing finite WKB curves). These curves will divide the 
punctured sphere up into cells with each cell denning a family of homotopically equivalent 
generic WKB curves as shown in figure 8 for an example of the 4-punctured sphere. To 
construct the triangulation, choose a representative curve from each cell, e.g. any one of the 
silver curves shown in figure 8. The claim is that the collection of these representative curves, 
which we will call edges, gives the desired triangulation [12]. 10 As a concrete example, the 
triangulation associated with the cell-construction of figure 8 is shown in figure 9. This same 
triangulation will play an important role in the 4-point function computation below. 

10 To see this in general consider a single zero of T as shown in figure 7. The zero is on the boundary of 
three cells. Choosing edges from the family of curves in each cell we see that they form a triangle. Thus the 
edges form a triangulation of the punctured sphere with each triangle containing a zero of T. 
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p 4 Pi 

Figure 9: The WKB triangulation of the 4-punctured sphere following from the WKB foliation 
shown in figure 8. Each edge E a b of the triangulation is a representative from one of the families 
of homotopotically equivalent lines in each cell of figure 8. This triangulation will be of central 
interest in the 4-point function computation. 



We have now finished the discussion of how to construct the WKB triangulation for a 
given T and Arg (£). Before moving on the the next section let us discuss one final point. In 
the following it will be useful to lift edges of the triangulation to the double cover E and to 
endow the lifted edges with an orientation. Recall that to = y/Tdw is a single valued form on 
S. Let dt be a tangent vector of the lifted edge E at a point on S. There are of course two 
possible orientations for d t . Note that by virtue of (32) we have e~ tl ^u • d t e M. We define 
the orientation of the lifted edge E by the condition e~ l ^u ■ d t > 0. Notice that each edge 
on the punctured sphere will lift to two edges - one on each sheet of E and that these two 
edges will have opposite relative orientation. Picking a particular orientation of some edge 
is equivalent to picking a particular sheet of E. 

4.5 Coordinates 

From the WKB triangulation we will now construct the so-called Fock-Goncharov coordi- 
nates [12]. These are natural objects to work with because they are gauge invariant and 
independent of the normalization of the small solutions. From the coordinates we will be 
able to extract the //-cycles that we need to compute the action (18). 

Consider some edge E of the triangulation. This edge is shared by precisely two triangles, 
and these triangles form the quadrilateral Qe (see figure 10). Number the vertices of Qe 
as shown in figure 10 with E going between P x and P 3 . As we mentioned in section 4.1, 
associated with each puncture P a there is a small solution s a . The solutions cannot be made 
globally smooth and single valued on the punctured sphere due to the monodromy around 
each puncture. However, we can define them such that they are single valued and smooth 
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Pi 

Figure 10: The two triangles sharing the edge £13. These two triangles define the quadrilateral 
Qei 3 , which is shown in gray shading. Each blue dot represents a puncture, which are the verticies 
of the triangulation and each black line and is an edge. 



throughout Qe- We then define the Fock-Goncharov coordinate as [12] 

, Asi A s 2 )(s 3 A Si ) 

Xe -[.-!)-, T7 r (odj 

(s 2 A s 3 )(s 4 A si) 

where all the s a are evaluated at a common point in Qe- 

As a concrete example consider the triangulation of the 4-punctured sphere shown 
in figure 9. In figure 11 we show how to apply the procedure just described to construct 
the coordinates corresponding to edges E 2i and E 2 ~ 4 . Consider first the left panel of 11. 
We define each solution s a throughout Q24 by parallel transporting from each P a where the 
explicit form of the solutions is known - see (21). The red lines indicate the parallel transport 
of each s a from P a to a common point A; clearly we can define the small solutions at any 
point A e Q24 in this way. Further, if the paths never leave the quadrilateral (or at least 
is always homotopically equivalent to a paths that never leave the quadrilateral) then the 
solution denned in this way is guaranteed to be single-valued and smooth throughout the 
quadrilateral, as required. With the solutions defined at a common point in the quadrilateral 
we can construct the coordinate X24, which is independent of the choice of A e Q 2i . Now 
consider the right panel of figure 11 where the grey shading indicates the quadrilateral 
associated to edge Q 2 \- These figures should be imagined on the sphere. Now to transport 
the small solutions to a common point one cannot avoid passing under a cut onto new sheets 
of some of the small solutions. For example s 2 must pass onto a new sheet in order to be 
smoothly continued to the point A. This is because if we were to compare the s 2 of the left 
panel and the s 2 of the right panel (by moving each respective A to a common point A 1 along 
the edge for example) the two paths of continuation would differ by a holonomy around 
P3, and thus the values at the point A 1 would not coincide but would differ by the action 

n We will show this in some concrete examples momentarily. 
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Figure 11: Here we show how to construct the coordinates X24 (left panel) and Xi 4 (right panel) of 
the triangulation of figure 9. The gray shaded regions represent Q24 and Q 2 ~ 4 respectively. These 
figures should be pictured on the sphere. The dashed blue lines emanating from the punctures 
indicate our conventions for defining the sheets of the small solutions as explained in section 4.2. 
The red lines indicate how we globally define the solutions s a by transporting away from P a using 
the connection. We use paths that never leave the quadrilateral such that the solutions used to 
form the coordinates are guaranteed to be single-valued and smooth throughout the quadrilateral, 
as required. 



of M 3 1 . Of course which solution we call S2 and M 3 ±1 S2 is purely a matter of convention. 
Similarly, which solutions acquire factors of M a depends on the choice of the point A. We 
stress that the coordinates are independent of all such ambiguities, as one can easily check 
using identities such as (M c s a A Sb) = (s a A M~ 1 Sb), etc. Then from figure 11 we read off 

, n (s 2 As 3 )(s 4 Asi) (M^ 1 S2AM 4 s 1 )(s 4 As 3 ) ( ^ 



S2 



We will also need to construct coordinates in the slightly degenerate case where the 
coordinate corresponds to an edge ending at a vertex that has only two incident edges 
(including the edge under consideration) for example all edges in figure 9 except E24 and 
E 24 . We show how to construct this coordinate in figure 12. Using the procedure described 
there we find 

(s 2 A M 3 s 4 )(s 3 A s 4 ) . (si A M 1 ~ 1 g 4 )(g 2 A s 4 ) 

X23 ~ [ } (M 3 s 4 A s 3 )(s 4 A s 2 y XU ~ [ \Mr 1 s 4 As 2 )(s 4 As 1 ) 1 J 

The other two coordinates X34 and X14 are computed in a similar way. 

We have now completed our discussion of how to construct the coordinates. Before we 
continue, let us comment on a useful property of these objects. Consider multiplying all of 
the coordinates associated with edges meeting a given puncture P. For example, the edges 
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Figure 12: Here we describe the construction of the coordinate for the slightly degenerate case 
where the coordinate corresponds to an edge ending at a vertex that has only two incident edges 
(e.g. P2 has only 2 incident edges: £12 and £23)- We construct the coordinate for edge £23 of the 
triangulation shown in figure 9. The quadrilateral prescription described above still applies, but 
one must take care to correctly define Qe and the solutions inside Qe- First of all, in order to have 
single- valued and smooth solutions throughout Q23 we must exclude a region between £3 and £4. 
Otherwise Q23 would contain £3 and thus the solutions could not be single valued in Q23 (there 
would be a monodromy around £3). Since the boundaries of the quadrilateral must be edges of the 
triangulation, the only choice is to remove a thin region running along edge £"34 and then to treat 
the two 'sides' of £34 as different edges. In the figure we have represented this process by showing 
£34 as doubled and with the region between the new edges excluded from Q23- We then define 
the solutions throughout Q23 in the same way as described in figure 11, by analytically continuing 
the solutions along paths from P a to A that stay within Q23 which is represented as the shaded 
region. Once we have defined the solutions at a common point we form the coordinate X23 given 
in equation (35). 

ending at P 2 in the triangulation of figure 9 are Ei 2 , E 2 ~ 4 , E23 and E24. Using (34)-(35) we 
have 



This property is true in general since the inner-products in the coordinates telescopically 
cancel in the product and the only thing that remains is the effect of the monodromy around 
the puncture which produces a \j? p factor. Thus we have the general rule [12] 



X12X24X23X24 = H 2 ■ 



(36) 



(37) 
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(s 2 A s 3 ) ~ i'3 2tJ /5 




Figure 13: Computing the £ — )• asymptotic of the coordinate X24 for a typical WKB triangulation. 
The blue disks represent the punctures and the black lines represent edges of the triangulation. A 
yellow x represents a zero of to and the wavy yellow line shows our convention for defining the 
sheets of E. The black arrows running along the edges indicate the choice of the direction for the 
edges. Each red curve indicates the transport of a small solutions in the limit £ — )• 0. The dashed 
red lines correspond to the transport of a solution appearing in the denominator of the coordinate. 
The transports used to form the coordinate combine into the continuous integral of to near the 
boundary of Q24; which can then be deformed into the cycle integral 724 shown in gray. 

4.6 WKB asymptotics of the coordinates 

The advantage of using the WKB triangulation is that the £ — > 0, 00 asymptotics of the 
coordinates of the triangulation are easily extracted given the discussion of section 4.3. 
That is, because we have maximum control over the large/small £ asymptotics of the small 
solutions when we transport along WKB curves. We give only the basic idea of the derivation 
of these asymptotics here and refer the reader to appendix B and [12] for details. 

To obtain the asymptotic of some xe one simply needs to use expression (30) for each 
inner-product of the coordinate, taking care to account for the direction of the WKB lines. 
Consider the coordinate associated with edge P24 in figure 13. The expression for this 
coordinate in terms of the small solutions is given in (34). We will now use formula (30) to 
compute the asymptotic of this coordinate in the £ — > limit. Let us take the directions of 
the WKB lines to be as given in figure 13. To evaluate the inner product (S2 A S3) we must 
transport the solutions to a common point. Since there is a WKB line flowing from P3 to 
P 2 we can safely use (to leading order) expression (30) to transport S3 to the neighborhood 
of P 2 , giving a contribution of the form (s 2 A s 3 ) ~ es Js To evaluate (s 3 A s 4 ) we must 
transport S3 to P4 since that is the direction of the WKB line and thus we get the contribution 
(S3 A S4) ~ ea^ w e may then reverse the order of integration in (53 A S4) and also move 
it to the numerator of the coordinate. Then the integrations from (S2 A S3) and (S3 A S4) 
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combine nicely into a continuous integral running just inside the boundary of Q24 from the 
neighborhood of P2 to P3 to P4. Repeating this analysis for the remaining two brackets one 
obtains a closed cycle integral passing along the boundary of Q2A- Recall from the discussion 
of section 4.4 that each triangle in the WKB triangulation encloses one zero of u. The 
integral of u thus encloses two zeros and so it can be deformed to the cycle integral shown 
in figure 13. Thus the non-vanishing contribution in the limit £ — > is given by 



Xe ~ (-l)exp -f 



u + Cf 



IE 



(38) 



The contour 7^ is the cycle encircling the two zeros contained in Qe and its direction is 
the same as that of the WKB lines corresponding to the brackets in the numerator of the 
coordinate. The term C$ is the O(£ ) contribution to the WKB expansion, which we will 
discuss momentarily. The overall (—1) prefactor in (38) is the same (—1) appearing in the 
definition of the coordinate (33). 

To derive the subleading WKB corrections (in the £ — > limit, for example) is essentially 
a matter of perturbation theory once the singular contribution has been extracted. We give 
a detailed discussion of this in appendix B. Here we will simply focus on the result and its 
implications. We find the first subleading contribution is given by 



C 



(0) 



log( 



± iir 



where ue is the number of w-spikes enclosed by 7^. 

Finally the £ — > 00 asymptotic follows in the same way as the £ 
integral around Qe of £w. 

To summarize, the £ — > 0, 00 asymptotics for \e are given by 



(39) 



and leads to a cycle 



X E~(-l) UE exp 



L J IE 



(40) 



where je is the cycle encircling the two zeros contained in Qe and its direction is the same 
as that of the WKB lines corresponding to the brackets in the numerator of the coordinate. 
Now it is clear how the choice of spikes (i.e. the choice of signs in (11)) is encoded into the 
coordinates - via the constant term in the WKB expansion which contributes the (— 1) UE 
factor in (40). Recall that ue is the number of w-spikes encircled by 7^. 



4.7 Shift relation. 

In section 4.1 we explained that there are two special solutions sp, Sp associated with each 
puncture P and that they are related to each other by a shift in the spectral parameter: 
sp(£) = a 3 sp(e~ t7T ^). Here we give an alternative relation between the small and big solutions 
that does not involve shifting the spectral parameter. The solutions sp and sp are linearly 
independent and thus we can expand any solution sq in terms of them. In particular we 
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have 

(s P A S Q \ ( S P A 8q\ _ 

SQ = {jp-^p) Sp+ {jp-A-i-p) Sp (41) 

M pSg = ( y^i) ,PS P + ( S -^) tf Sp (42) 
\SpAspJ \spAspJ 

For the second equality we have used (27)-(28). Combining these two equations it follows 
that 

' M PSQ As Q \ /gpA^\ 

MpSqAsp) y PP '\spAsp) v ; 

The utility of this equation is that it allows us to replace certain wronskians involving big 
solutions (as on the RHS of (43)) in terms of small solutions with monodromies. This will 
play a key role in the derivation of the functional equations that we present in the following 
section. 



4.8 x _s y s tem. 

We will now derive a set of functional equations for the coordinates which, together with cer- 
tain analytic properties, allows us to determine the coordinates completely. Our inspiration 
comes from the solution of the bosonic Wilson-loop problem at strong coupling [10] where 
the solution involves a set of functional equations of the schematic form 12 

Y+Y- = F a {Y) (44) 

where f ny± = /(#±m7r/2). On the RHS of (44) the function F a can depend on all of the Y a , 
but with their arguments un-shifted. The only shifts in the spectral parameter occur on the 
LHS of (44). For the Wilson-loop problem the F a are such that (44) takes the form of a so- 
called Y-system which commonly appear in the context of 1 + 1 dimensional integrable QFT's. 
Here, using the general formalism of [12], we will arrive at a set of functional equations with 
the same schematic form as (44) but with the F a of a different form than that occurring in 
the Wilson-loop problem. We will call this type of functional equation a x-system. 
To derive a relation of the form (44) we begin with the LHS. Using (24) we have 

XeXe = XeXe + ( 45 ) 

where \e is defined by taking \e and replacing each small solutions s a — » s a . To obtain 
the schematic form (44) we need to rewrite (45) in terms of only un-shifted small solutions. 
That is, we need to get rid of all the tildes without introducing any shifts in the spectral 
parameter. For this we can use (43) after applying the Schouten identity 13 to (45) to obtain 

XEX E {1 + Abc )(l + A da ) (46) 



12 The linear problem associated with that problem is very similar to the one considered here and the Y a 
are (up to shifts in the spectral parameter) the coordinates associated with that problem. We are referring 
here to the special case where the Wilson loop lives in an IRi.i subspace. 

13 (s a A s b )(s c A s d ) + (s a A s c )(s d A s b ) + (s a A s d )(s b A s c ) = 0. 
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Figure 14: Graphical rules for constructing £ (P; Q — > Q). Start at edge Ppq and continue in a 
counterclockwise fashion about P forming the nested product (49) by multiplying the coordinates 
for each edge encountered along the way (i.e. the coordinates associated with each edge intersected 
by the red line in the order indicated by the arrow). The dashed blue line indicates our convention 
for cutting the solutions to account for the monodromy around P. The small solutions used to 
form the coordinates are defined in the vicinity of P by analytically continuing them throughout 
the triangles along the direction indicated by the red arrow and thus if we use sq in XP,a then we 
must include a monodromy matrix when the solution is continued around P to form XPl- 



where we have defined the useful auxiliary quantity 

(sq A S P ) (s P A Sq 
(Sp A Sp) (s Q A Sq 



a P q = (-i) { ; qA : p I { ; pA I q \ m 



-1) (1 + (1 + ( MpSq*sq \ ( M Q spAs P \ 

^ P> K PQ> \MpSqASpJ\MqSpASqJ v 



Here, the edge E is the edge ac in Qe where the vertices are labeled abed in counter-clockwise 
order. To go from (47) to (48) we used the shift relation (43). The last step is to rewrite 
the wronskians appearing in (48) in terms of the coordinates. Once this is done, combining 
(45) — (48), we can assemble a functional equation of the form (44). To do this (following 
[12]) we introduce the quantity 

£ (P; Q -+ Q) = 1 + X P,a (1 + XP,a-i (1 + ...Xp,2 (1 + Xp,i))) (49) 

The coordinates appearing in this object are shown in figure 14. By repeatedly applying the 
Schouten identity (see footnote 13) starting with (1 + Xp,i) one can see that the Wronskians 
in (49) telescopically cancel so that 14 

yj /p. q y q\ __ (go A gq+i)(sp A s a ) = (s P A s a )(Mps Q A s Q ) 
(so+i A s a )(s A s P ) (sQAs a )(Mps Q As P ) 



14 An easy way to see this in general is to use induction [12]. The case a = 1 is simple to prove using Schouten 
identity. Then one can show (again using Schouten) that £ (P; Q a +2 ~ > Qo) = 1 + Xp,a+iS (P; Q a +i Qo)- 
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In going from the first equality to the second in (50) we have accounted for the monodromy 
acquired by the small solutions when they are analytically continued around P (see figure 
14). Then, from (50) and (48) we have 



Finally, using (51) in (46) and noting (49) we obtain a closed functional equation for xe of 
the form (44). Repeating this procedure for the coordinate associated to each edge in a given 
triangulation gives the desired set of functional equations. Note that this procedure can be 
applied to derive the x-system for an arbitrary number of punctures. In section 5 we will 
apply this procedure to the triangulation (9), which is one of the triangulations of interest 
for the four-point function computation. 

4.9 Inverting x-systems 

In the previous section we showed how to derive the x-system associated with a given trian- 
gulation of the iV-punctured sphere. In this section we will discuss how to use the x-system 
along with certain analytic properties of the coordinates to obtain integral equations that 
determine the xe uniquely. 

The basic idea behind the inversion of a x-system is to Fourier transform (the log of) 
each equation since in Fourier space these nonlocal relations become local as the shifts in 
the parameter 6 can be undone in the usual way. For such a procedure to be successful one 
must have a certain amount of control of the analytic properties of the coordinates. Let us 
discuss this carefully. The equations that we want to Fourier transform have the form 



where Fe (x) is an explicit function of the coordinates which follows from the discussion 
of section 4.8. We have introduced the arbitrary shift <p for reasons that will be explained 
momentarily. Note that XeXe + = XeX~e~ since the small solutions are 27ri-periodic, which is 
why we can have either shift Fe (x ) on the RHS. The choice of this shift is arbitrary since 
the objects we will eventually compute (the ^-cycles) are functionals of the coordinates only 
through Apq which is «7r-periodic and thus does not care about the choice of shift. As a 
convention we choose the shift —m/2. 

To Fourier transform the relationship (52) one must be sure that the transform converges. 
Moreover, to undo the shifts on the LHS, one must account for the singularities (if any) of 
log xe in the strip of width it centered along the line where the transform has been performed. 
We will now discuss each of these issues in turn. 

The information from the WKB analysis will allow us to ensure the convergence of the 
Fourier transform, provided certain conditions are satisfied. First consider the LHS of (52). 
We need to ensure that the transform of each individual term converges. We can ensure this 
if we know the asymptotics of the coordinates in the full strip Im (9) G (0 — 7r/2, <p + tt/2). 
The coordinates should be derived from the triangulation that one has at Im (6) = 0. Then 
the WKB analysis guarantees that the asymptotics are given by (40) in a strip that includes 
the region Im e (<f) — 7r/2, <fi + 7r/2). Each term on the LHS can be made safe to transform 



(1 + fi P ) 2 (l + fx Q ) 2 A PQ = xpqS (P; Q^Q)E(Q;P^P) 



(51) 



log X E (0 + i( f>) + lo S X E ( d + i( f>) = lo S f e (x ± (0 + #)) 



(52) 
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by making (on the LHS only) the replacement xe — > Xe/x^e where Xe * s the asymptotic 
(40). This replacement does not modify the equation since \ Xej {^ej = ^' 
Now consider the RHS of (52), which has the form (see equation (46)) 



logF E ( x ± (9 + i(f>)) =log 



:i + A abK i +Acd ) , n±iK/2 ^ 



;i + A bc )(l + A da ) 



(53) 



Each Apq is computed by (51) and (49). For the RHS of (53) to be decaying it is sufficient 
for all of the Apq in (53) to be decaying. If all the x-functions are decaying then from (51) 
and (49) it is clear that all of the A P q will decay; the //-factors will decay by virtue of the 
rule (37). On the other hand, if all the x-f unc ti° ns are growing the //-factors in (51) will 
dominate the RHS of (51) so that Apq is still decaying; to see this one should re-express 
the //-factors in terms of the coordinates using (37). Thus the RHS of (53) will decay if 
all of the ^-functions are growing, or alternatively if they are all decaying. For generic it 
will generally not be true that the RHS of (53) is well behaved, and one must try to find a 
range of 0-values for which the xe are all decaying or are all growing. If a suitable can be 
found, then (52) can be directly solved by Fourier-transform. In all of the examples we have 
considered (in particular, those relevant for the 4-point function) it has been possible to find 
such a <p. 

Concerning the issue of singularities within the strip of inversion, it follows from (19) 
that the Wronskians (s a A s&) (9) are (in an appropriate normalization) analytic away from 
9 = ±oo. It is, however, possible for these objects to have zeros and in the following it is an 
assumption that there are no zeros in the strip where we do the inversion. 15 In section 5.1.4 
we perform numerical tests that support this assumption. 

Finally, we use the Fourier analysis to obtain 

where Xe{9) = Xe (9 + i<p — iiv/2) and Xg is the (shifted) asymptotic (40) and Fe(X) 
is an explicit function of the coordinates which follows from the discussion of the previous 
section. 

The equations (54) can easily be solved for the Xe by iterating them in a computer. In 
the next subsection we will show how to extract the //-cycles of formula (18) from the Xe 
which are computed using (54). We will then perform some numerical tests in section 5.1.4. 



4.10 Extracting 77-cycles 

Once the coordinates are computed according to the prescription of the preceding section 
we extract the r/-cycles as follows. What we need to compute are the individual Wronskians 

lD In the limit where the WKB approximation holds, i.e. when — > ±00 or in the limit of large zero 
modes \Ze\ — > 00 [12], it is clear that (in an appropriate normalization) the Wronskians will not have any 
zeros since (suppose we compute the Wronskian near Pb) then s a will be the big solution near Pb and is thus 
linearly independent of s& which is small at Pb. For finite values of 9 (or alternately of \Ze\) we have no 
concrete way of arguing that these zeros are not present. 
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(s a A s b ). For this, note that from (47) and footnote 13 we have 



:i + a 



ab ) 



[S a A s b ) (s a A s b ) 



(55) 



We can choose a guage where (sp A Sp) 
With this gauge choice we have 



(S a A S a ) (Sb A S&) 

1. The final result will be gauge independent. 



log (s a A s b ) + log (s a A s b ) + = log (1 + A ab ) 



(56) 



Here we will use the notation + 20 where 9 and <fi are real. We then insert the 

zero-modes on the LHS in the same way as for the x-system (see section 4.9). We are only 
interested in P a and P b that are connected by a WKB line when Arg(£) = 0, and thus 
we have good control over the asymptotics in the required strip. Performing the Fourier 
transforms we obtain 



log (s a A s b ) (6 + i(j)) = 

where we have defined 

w ab = lim lim 



b uj ab + ^e e+ ^w ab ) 



! Vfdw + ~y log(w a 



d9' \og(l + A- b (9' + t 
2tt cosh (6 - 6') 



(57) 



(5? 



The integration in (58) is performed along edge E ab . The direction of integration is the same 
as the direction of the edge E ab (see appendix B). Note that the logarithmic terms precisely 
cancel the divergence from the endpoints of integration in (58) so that the w ab are finite. In 
going from (56) to (57) we have used the asymptotics for (s a A s b ) derived in appendix B. 
Expanding (57) around 9 — > — oo, and comparing with (26) with £ = e e+l< ^ we read off 



(W 
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log(l + A a " fe (^ + z0)) 



(59) 



The contour of integration in J E rj is along the WKB line connecting P a and P b and the 
direction of integration is the same as the direction of the edge E ab . This formula allows 
us to compute the //-cycles from the x-functions since the Apq are explicit functions of the 
coordinates. 



5 The AdS action 

5.1 Regularized AdS action 

Now that we have introduced the needed tools we are ready to calculate the action (18). We 
will demonstrate for the case of the 4-point function, but the method is general and could 
be performed for any number of operators inserted along a line. The computation will be as 
follows. First we will introduce the relevant WKB triangulation which will be topologically 
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equivalent to the triangulation shown in figure 9. Second, using the procedure of section 4.8 
we will derive the x-system satisfied by the coordinates of this triangulation. Supplementing 
these functional relations by the WKB asymptotics we will invert these functional relations 
using the technique of section 4.9 to obtain a set of integral equations that uniquely determine 
the coordinates. Finally, from coordinates we extract the 77-cycles using the method of section 
4.10. Once we have the //-cycles, we compute the action using (18). 

5.1.1 Stress-energy tensor and WKB triangulation 

For the purpose of the following computation, a useful parameterization of the stress energy 
tensor is 

_ 1 / c + CiW + c 2 w 2 + Uw 3 
[W) ~ (w - w 4 ) 2 \ C °° (1 + wf{\ - wf 

Here we have fixed three of the insertion points at W\ = +1, W2 = 00, W3 = — 1 using 
the world-sheet conformal symmetry. The fourth insertion point is left at the position w 4 
which should be fixed at the saddle point w\ once the full action is assembled. For the 
purpose of demonstration we will take W4 to be between u> 3 = — 1 and W\ = +1. When the 
dominant saddle point is located in one of the other intervals one can proceed by a similar 
procedure. The constants c a = c a (w&, A) are functions of W4 and dimensions of the operators 
and are fixed by the condition (7). Their explicit expressions are given in appendix D. The 
parameter U is unfixed by the condition (7) and implicitly parameterizes the cross ratio of 
the four operators (recall that they are inserted along a line in the boundary theory so that 
there is only one cross ratio). The analytic structure of T, the resulting WKB-structure and 
the WKB triangulation are shown in figure 15. 




5.1.2 ^-system for the 4-point function 

From equation (46) and figure 15 we have 

X24X24 - { X ^J - (l + A 34 )(l + A 12 ) (61) 



To compute each Apg we use formulas (51) and (49) along with the rules given in figure 14. 
In that way we find 

A X24 (1 + Xl2 (1 + Xi4 (! + X23))) (1 + X43 (1 + XA2 (! + X4l))) 

Mi - (1-^)(1-^) (63) 

A X23 (1 + X34) (1 + X24 (1 + Xl2 (1 + X24))) If. AS 

= MIR) (64) 

with the rest of the Apq being related by relabeling (see appendix D for the explicit formulas). 
These expressions and equations (61) — (62) provide a closed system of functional equations 
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Figure 15: Constructing the triangulation for the 4-point function. In the left panel we show the 
WKB cells for Arg (£) = 0. The cell walls are formed by the separating WKB curves as described 
in section 4.4; as described there, inside each cell there is a 1-parameter family of generic WKB 
curves and by taking a representative curve from each family we obtain the triangulation shown in 
the right panel. In the right panel the black lines are the edges of the WKB triangulation and the 
wavy yellow lines show our convention for defining the branches of to. Notice that this triangulation 
is topologically equivalent to the one shown in figure 9. This means that we can borrow the results 
derived for that example. In particular, the coordinates can be carried over from that example by 
making the proper identifications. The cycles corresponding to each coordinate are represented by 
the gray curves - we show only the portion of each cycle on the sheet of oj where the edge -E34 has 
orientation towards P4 as indicated by the black arrow along edge -E34. 



for the 6 coordinates associated with the triangulation shown in figure 15. 

These functional equations can be converted into integral equations of the form (54) using 
the technique described in section 4.9. To apply the procedure of section 4.9 one must find a 
(f> such that the RHS of (52) is decaying, and for this one should appeal to the WKB analysis. 
The WKB cycles which determine the asymptotics of the coordinates are shown in figure 15. 
When Ai ~ A 3 and U ~ 0, ~ the cycles shown in figure 15 all have Arg(_^> u) ~ 7r/2. 16 
In this case = is a suitable choice since then all \~e w * n be growing and (53) will decay 
rapidly. 17 In summary, the integral equations in the region of present interest are given 
by equations (54) with Fe given by (61) — (64). These equations will remain valid for all 
values of the parameters A a , U, and W4 such that the triangulation is unchanged. If one 
deforms these parameters too much the triangulation will jump. One can then easily write 
the x-system for the new triangulation and apply the same procedure to obtain the integral 

16 Interestingly, when Ai = A3 and U — W4 = there is a symmetry which causes the RHS of the ^-system 
to trivialize (i.e. to reduce to 1 for all \e) and the ^-functions can be computed explicitly (they are just 
equal to their zero-mode part). This is reminiscent of the case for the three-point function and, in fact, there 
is also a change of coordinates that maps the specific case Aj = A 3 and U = W4 = to two copies of a 
three-point function. 

17 This will continue to be the case as long as the Arg(J> uS) remain in the upper-half plane. In other- 
words, the inversion procedure will be valid for all U and W4 such that the triangulation is unchanged since 
the triangulation will jump precisely when one of the § u) crosses the real-axis [12]. 
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equations for that region of parameters. 18 

By numerically iterating these equations (using \e as t ne initial iterate for each \e) 
we obtain the x-f unc ti° ns - The //-cycles are then extracted from the x-f unc ti° ns using the 
procedure of section (4.10). In the following section we will write the regularized AdS action 
in terms of these //-cycles. 



5.1.3 Finite part of AdS action 

Now that we are able to compute the //-cycles (see previous subsection) we can use the 
formula 

A fm = jf Vrf (cosh 7 - 1) = I - \ (Ju^ I-, 1 (jf 77) . (65) 

(see section 3.2 and equation (18)) to compute the regularized part of the AdS action. To 
use (65) there are few steps. These steps are simple but tedious and we will only list them 
here (see appendix D for a detailed implementation). As described in section 3.2 one should 
first modify T by spreading the double poles slightly such that u = y/Tdw has an additional 
square-root cut at each of these points. Then one should choose a complete basis of a- and 
6-cycles (five of each is needed for the 4-point function). One can then apply formula (18) 
and then take the limit in which the small cuts close to form simple poles in u>. Once this is 
done the area will generically be expressed in terms of three different types of //-cycles: cycles 
connecting two punctures, cycles connecting a puncture with a zero and cycles connecting 
two zeros. The latter two can be expressed as linear combinations of the puncture-puncture 
cycles as described in appendix D. Once this is done, the final result takes the elegant form 

A fin = I ~i y^UEVE (66) 
EeT 

where the sum runs over the edges in the triangulation (see figure 15), rjE ab is defined in (59) 
while UE ab is the w-cycles that intersects edge E a b (i.e. the integral of u that is associated 
with the coordinate Xab] these integrals are shown as the gray contours in figure 15). 19 

Formula (66) and the procedure of section 4 for computing the recycles solve the problem 
of computing the regularized AdS contribution to the 4-point function. In the next section 
we present some numerical tests of the procedure. Let us note that the procedure of section 
4 is general and can be implemented for any number of punctures. Further, while we have 
only proved equation (66) for the case of the 4-point function, given its simplicity one might 
suspect that the formula holds in general (with tt/3 — > n/12 x (^number of zeros of T), of 
course). 20 Even if the general result does not take the simple form (66), for a given T (i.e. 
for any number of punctures) the procedure described in section 3.2 is still valid and one 
can still write Afi n in terms of the t\e for the corresponding triangulation). In principle this 

18 Another (more elegant) approach would be to find a systematic way of analytically continuing the integral 
equations from one region of parameters to another as was done for the TBA equations of [10]. 

19 Note that in formula (66) both integrals uje and tje are the segment integrals between the appropriate 
limits. For example, the cue = \ § w. In this sense we are abusive with the term 'cycle'. 

20 It would be a simple matter to check this, but we have not pursued this issue. We did check that the 
formula holds for the 3-point function (see appendix E). 
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A 3 


A 4 


Ai 


A 2 


Numerics 


X-system 


1/5 


1 


2 


1 


2 


0.84807 


0.84812 


1/2 


1 


2 


1 


2 


0.82421 


0.82423 



Table 1: Comparison of the Afi n obtained by numerically integrating (10) and the area computed 
from the x-system. The results are for the spike configuration of figure 4B. 



solves the problem of computing the regularized AdS contribution to the iV-point function. 
We have performed numerical tests only for the case of the 4-point function. We present 
these numerical results in the following section. 



5.1.4 Numerical tests 

We now present numerical tests of the method described above. We solved numerically the 
modified sinh-Gordon equation (10) for the function 7 and then using this numerical solution 
to directly compute Af in via 

J Vrf (cosh 7-I) (67) 

The general set-up of the numerical problem essentially follows that of [8] with some modifi- 
cations. However, the numerical method that we use to solve the PDE (10) is quite different 
from that of [8]. 21 We place the punctures at w 3 = —1, w 4 = 0, W\ — 1, and w 2 = 00. We 
then map the sphere to a square domain with the point at infinity mapping to the boundary 
of the square and the real axis mapping onto itself. Since 7 must vanish at the punctures, we 
should impose 7 = along the boundary of the square domain since W2 maps to the bound- 
ary of the square in the new coordinates. Further, since for either configuration of spikes (see 
section 3.3 and figure 4) there is a fold-line along the real axis, we know 7 (x, 0) = where we 
are using the coordinates w = x + iy and writing 7 = 7(2, y). Thus we can solve the problem 
in half of the square with the Dirichlet boundary conditions 7 = on the boundaries. Lastly, 
we must remove the logarithmic singularities (11) in order to have a nice smooth function 
to solve for. A suitable function is 



27re 9 = 7 + 2 ° a log 



(lU - Zg) (W ~ 

(1 + ww) 



(68) 



where we a a = ±1 is determined by 7 ~ — a a \ log TT at z a . The numerator of (68) removes 
the log divergences (11) in 7 while the denominator is included to kill off these additional log 
terms at infinity. In the numerical implementation we fix the spike configuration we want 
to describe by choosing the set of {cr a }. Finally, to numerically integrate the equation (10) 
(re- written in terms of 7 re9 , of course) we use a standard relaxation method with an uniform 
grid. 

In table 1 we compare the numerical results with the analytic results. The numerical 



21 We are very grateful to Romuald Janik for providing us with a copy of the code used in [8] which was 
very useful in helping us to develop and test our own numerics. 
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Figure 16: Here we show the values of 7/14 and r/34 evaluated along several different contours. For 
example, the column labeled r/14 [A] ( [B] ) shows the values of 7714 for the spike configuration of figure 
■iA(B) for each of the contours shown to the right of the column. We use the parameter values 
A3 = Ai = 1, A2 = A4 = 2 and U = 1/5 for both spike configurations. There are five digits that 
we trust since they are unchanged for the different contours and they should be compared with our 
result from the functional equations that is rjf^ S ^ S em [A] —0.033169, r]^ sys ^ em [A] « —0.014503 

and r]X~ system [B} » -0.031628, T& aysbem [B] « -0.013588. In the digits where the forms are 
closed there is perfect agreement with the analytic results. 



results are obtained by the area computed using (67) with the numerical solution for 7. 
The analytic result is obtained from (66) with the 77-cycles computed using the x-system 
procedure. These results show a good agreement of our formula with the numerics. 

A sharper measure of the agreement between the analytics and numerics is to compare 
directly the 77-cycles. In figure 16 we show the numerical results for r/ 14 and 7734 computed 
along several different contours. This allows us to test the closure of the numerical rj which 
we obtain from the numerical 7 via (16). Note that closure of rj implies that 7 must obey 
(10) and thus this is a good measure of the numerical error. Indeed, one can see in figure 16 
that the numerical cycles agree with the analytical predictions in all digits in which they are 
closed. That is, the numerics is in agreement with the analytics in all of the digits for which 
the numerics can be trusted. 

Finally, it would be interesting to perform numerical tests for a larger portion of the 
parameter space (i.e. more values of the A , U and W4). To perform a systematic study 
will probably require an improvement of our numerical method as our current method, while 
extremely simple, has very slow convergence. 
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5.2 Divergent part 



In section 5.1 we completed the task of computing the first term in formula (14). In this 
section we will discuss the second term 

-/ d 2 wVrf=-^y2A 2 a \oge a -A reg (69) 

•M{e tt } 2 a 

where A reg is finite at e a — > 0. The contribution A reg can be computed by simple but tedious 
application of the Riemann bilinear identity and there are many ways to write the result. 
For example 

A reg = i 22 W E^E - « 2^24 ~ ^i 4 ){^2A ~ ^u) (70) 
EeT 

where the sum is over the triangulation shown in figure 15 and ^E ab = ^ab is defined in (58). 
The ue are defined in the same way as in (66). One can check this formula by comparing 
with the direct 2D numerical integration of VTT with small circular disks cut out around 
the puncture (in Mathematica one can use NIntegrate along with the Boole command, for 
example). 

We recall that (69) came from the regularization of the string action where we have added 
and subtracted VTT from the integrand of the AdS action. This integral depends explicitly 
on the cut-off e a around the punctures. It will be important to understand the connection 
with the physical cut-off £ at the boundary of AdS. Fortunately we can extract the needed 
information from the linear problem since we have good analytic control over the solutions 
near the insertion points. To proceed by this route (which parallels the discussion of [8] 
for the 3-point function) we must first describe how the string embedding coordinates are 
recovered from the linear problem formalism, which is via the aptly-named reconstruction 
formulas. We will discuss this in the next subsection, 5.2.1. After that, in section 5.2.2 we 
will use the reconstruction formulas to eliminate the e a in favor of £. From this procedure 
we will recover the standard spacetime dependence in (1) along with a contribution to the 
function f(u,v). This will complete the computation of the semiclassical AdS contribution 
to(l). 

5.2.1 Reconstruction formulas 

The reconstruction formulas allow us to express the string embedding coordinates in terms 
of solutions of the linear problem. This point is crucial in our construction for the following 
reasons. First, we have introduced some regulators in the world-sheet, e a , that must be 
related to the physical cut-off in the boundary of AdS, z — £. Second, by using them we will 
be able make the spatial dependence explicit in the final result, namely the insertion points 
x a of the operators in the gauge theory. 

Consider two solutions of the linear problem, ipA and tpB normalized as (ipA A^b) = 1, 
and construct a matrix \l/ as 

* = &a 1>b) • (71) 
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The matrix \1/ obeys the same equations of motion as ipA.B (19), namely 

(d + J w )^ = 0, (d + J^ = 0. (72) 
where J w and Jw are defined in (19)- (20). One can verify using (9) that the quantity 

?/ = ~tr (aVW*)| 0=Q (73) 

with a 1 = a 1 , a 2 = —ia 2 , a 3 = a 3 , satisifes the same equations of motion as Y 1 and also 
the constraint y ■ y = — 1 (with the AdS metric). In this way we establish a correspondence 
between target space coordinates and solutions of the linear problem, 

I = Y 2 -Y 1 = 2i* 1 i* 21 , - = Y 3 = i(q 11 ^ 22 + q 12 q 21 ) (74) 

z z 

In order to relate the operator insertion points x a and physical cut-off S with the linear 
problem data, it is convenient to express ipA and "0s i n terms of the elementary solutions s a 
and s a whose behavior close to the punctures is given by (21), 

IpA = (i>A A S ) S a + (s a A ^a) §a, IpB = {ipB A S a ) S a + (s a A l/) B ) §a (75) 

Close to the punctures the solution s a becomes dominant. Then, using (74) and the explicit 
form of s a close to the puncture P a we get that 

z = — ^— -^\w-w a \ Aa (76) 

% [s a A i) A ) 

where the subscript indicates that the solutions are evaluated at 9 = (recall that this is 
the value where the physical problem is recovered - see equation (73)). Equation (76) is the 
relation needed to make the connection between the world-sheet and physical cut-off's 

A a log e a = log £ + log | (s a A ipA) lo (77) 

Finally, using once again (74) we express the insertion points x a of the operators in the gauge 
theory as 

= 

5.2.2 Physical regulator and spacetime dependence 

We can now use (77) to eliminate the e a in (69) in favor of the physical cut-off at the boundary 
of AdS z = S. We have 



J2 K e a = I A a log S + A a log I (S a A IpA) lo j 
a \ a a / 



(79) 



where a and A refer respectively to the small solution s a and one generic solution ipA appear- 
ing in the reconstruction formulas. Now we will eliminate the factors | (s a A i^a) |o i n terms 
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of objects that we can compute. 

The terms | (s a A i/ja) lo can be related to the insertion points x a in target space and 
overlaps of the elementary solutions evaluated at 6 = through expression (78). Using 
Schouten's identity one can verify that 





(s b A s a ) 







As„) 





■^ba^ca 




A s b ) 






Sa A ip A ) lo = FT — 7 — v-j (80) 



for a, b, c distinct. This solution is unique up to different ways of rewriting the spatial 
dependence using the cross-ratio 

XU%23 (si A s 4 ) (s 2 A s 3 ) 

U = = r r— lot) 

X12X34 [si A s 2 ) {s 3 A s 4 ) 

where we have used (78). Note that we can compute the brackets appearing in (77)-(78) 
using (57). In particular we have 

A 1 \ f de log (1 + az) , s 

log (s a A s b ) = -w ab + -w ab + / 6 V , n abJ (82) 



2 au 2 "7 7 M 27r coshS 

This formula is valid when there is a WKB line connecting P a and P b . If a bracket appears 
for which we do not have a WKB line, we can simply use the cross ratio (81) to eliminate it 
in terms of brackets that can be computed using (82). 

Finally, using (80) in (79) and massaging the resulting spacetime dependence by extract- 
ing multiples of u and (1 + u) we find 

1 



= n(kA a6 |o)-^(^) ^ (83) 

a>b 

where A ab = (^2 A c ) /3 — A Q — A5. The extra factor of 2 in the exponent on the left hand 
side of (83) anticipates the sphere regularization which turns out to be similar to the AdS 
part and will be treated in section 6.1. 

We recognize in (83) the canonical spacetime dependence in the 4-point function of a 
conformal field theory (compare with equation (1)). The appearance of the cut-off in (83) is 
related to the renormalization of the operators. In fact, if we define 0& a = £ Aa O& a this will 
cancel the E factors in (83). To be more precise, we should define a 4-point function that 
is independent of the operator renormalization. The standard procedure is to divide by the 
appropriate product of 2-point functions such that normalization factors cancel. The same 
factors of £ will appear in these 2-point functions and will cancel with those in (83). We will 
thus drop the factors of £ in the formulas below. 



5.3 Summary of the AdS and divergent contributions 

We have now computed all the parts of (14). In this section we summarize the full result. 
The semiclassical limit of the 4-point function (1) is given by 

(fftff^SfkYilM**, (84) 

a<b 
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where the * denotes evaluation at W4 = w\ and we define 

//If M = e-^ A f- (85) 

/ d 1f XS M = e- 2 ^-n(|, a A Sfe | )^ A - (86) 

a>b 

and /j in will be defined momentarily. The contribution Af in is given by (66), A reg is given 
in (70), the brackets in f^ SxS are given by (82). 

The sphere part of the correlation function cont ains divergences of the same type as 
AdS. We therefore regularize it also by subtracting VTT. Such finite contribution is what 
we denote by fj in 

= g-^ /e( S5 contribution- VTt) _ ^ 

where S* 5 contribution stands for the 5* 5 Lagrangian and wavefunctions [8] . To compensate 
this subtraction, we include the factor of 2 in front of A reg in expression (86). In general we 
cannot complete the construction of the 4-point function because we are unable to compute 
the contribution ff in . Fortunately, for correlators involving only BPS operators of the same 
type (e.g. only Z and Z) the sphere part is known and we can assemble the full result. This 
is the subject of the next section. 



6 Full correlation function for BMN operators 

In this section we compute the full correlation function for operators of the type tr Z A when 
A scales as a/A. For these type of operators, the sphere part ff in was already known [31] and 
therefore we can complete our computation. We stress that, unlike the three point function, 
this four point correlator is not protected. In section 6.2, we fix the location of the puncture 
w 4 by the saddle point method and discuss some issues on the multiple string embedding 
configurations. In section 6.4 we perform an analytical check of our procedure by studying 
the extremal limit where A 2 = A 4 + A 3 + A 4 , which is known to be protected from quantum 
corrections. 



6.1 Sphere part 

The sphere part of the correlation function involves the classical wave-functions associated to 
the external states. We consider specifically the correlation function of four BMN operators 22 

(tr Z Al ( Xl ) tr Z k \x 2 ) tr Z^ 3 (x 3 ) tr Z Ki (x 4 ) ) , (88) 

for which the wave-functions are known [13, 14]. The string dual of these operators cor- 
responds geometrically to a string that is point-like in the sphere and rotates around an 
equator [15]. The surface developed by the worldsheet is not extended in the sphere. 

22 We are using the following notation for the dimensions of the operators A a = V\A a . 
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Let Xi (i = 1, 6) be the coordinates in S 5 . This particular string state can be expressed 

as 



Xt + iX 2 = e l 



Xi = 0, 



3,. ..,6 



(89) 



where p is an azimuthal angle of the sphere. The wave-functions for tr Z Aa and tr Z Aa are 
given respectively by 



JA a <p(w a ,Wa) 



-iA a ip(w a ,Wa) 



(90) 



where the field ip is evaluated at the puncture corresponding to the respective operator in- 
sertion. 

As the wave-functions scale exponentially with \/X, they will act as sources for the equa- 
tions of motion for p. The total sphere contribution is then given by 



exp 



7T 



(^j d 2 wdpdip + in (A 3 <p w= _x + A 4 p w=W4 - Ai^ w= i - A 2 (p w=00 )j 



(91) 



Considering both the contributions from the S 5 action and wave-functions as an effective 
action, we obtain the equations of motion for ip which are solved by 



<p(w, w) = i (A 3 log \ w + 1| + A 4 log \ w — w^l — Ai log \ w — 1|) 



(92) 



This solution has an additional singularity at infinity with charge —A3 — A4 + Ai(= — A2), 
corresponding to the wave-function inserted at infinity. This is consistent with i?-charge 
conservation. We may now plug (92) into (91), introducing cut-off's around the punctures 
to regulate this contribution. This amounts to evaluate the solution at a distance e away 
from the punctures. As in the case of the AciS action, the logarithmic divergences 



exp 



£A?lo. 



(93) 



need to be regularized. We do this by subtracting vTT from the integrand. To compensate, 
we add a similar contribution to the divergent part, that was already treated in the previous 
section (indeed, this regularization procedure is responsible for the factor of 2 appearing 
in front of A reg in expression for f^ SxS , see (86)). The dependence on the cut-off's then 
disappears yielding the following expression for the regularized sphere action and wave- 
functions 

K-i| AiA4 ^ 



./ 



s 

fin 



exp 



AU reg -log2 A3Al -log 



W 4 + 1| A 3^4 



(94) 



where ff in was defined in (87) 



6.2 Saddle point determination 

We have shown how to compute the quantities (85)- (86) as a general function of W4. However, 
to compute (84) we must evaluate at the saddle point W4 = wl- Before discussing how to 
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(a) 



(b) 



Figure 17: In (a) we computed several points for a region around u ~ 0.9 and interpolated the 
surface in blue. We then took a constant cross ratio slice. The intersection is given by the orange 
line. In (b) we have plotted the area as function of for cross ratio u = 0.9, where we have a 
saddle point at w\ « —0.24. 



locate the saddle point for generic values of the parameters, let us discuss a very symmetrical 
situation in which we can guess its position. 

Using conformal symmetry, we can fix three of the points in the target space at £3 = 
— 1, X\ = 1 and x 2 = 00 and also in the world-sheet at w 3 = — 1, W\ = 1 and w 2 = 00 . The 
fourth point 24 will then be related to the cross ratio. For the particular choice of the cross 
ratio (81) equal to 1, the fourth point will be located at zero. Moreover if we choose the 
points at X3 = —1 and x\ = +1 to have the same conformal dimension and the same type of 
fields (say Z's) then this is a very symmetrical configuration. Going back to the world-sheet 
coordinates, by symmetry we expect to find the saddle point at w* 4 = 0. 

Let us now describe how we determine the saddle point in practice. All the physical 
information that we input is contained on the stress-energy tensor. Besides the conformal 
dimensions of the operators and the position of the punctures, there is the extra parameter 
U that translates the additional degree of freedom of the cross ratio. The two are implicitly 
connected by the expression 



by formulas (81) and (34). Therefore, our strategy will be to compute the area and cross 
ratio for many points in a region in the (7/74, U) plane and then translate that information 
to the (w4,u) plane. Finally, we take slices of constant cross ratio, and determine for which 
W4 the area is stationary. 23 We can very easily perturb away from this very symmetrical 
configuration and track the location of the saddle point. For example, in figure 6.2, we have 
found a saddle point for u = 0.9 at wl ~ —0.24. One can compute the location of w* 4 to 
arbitrary numerical accuracy by iterating the x-system; here we present only 2 digits since 
we just want to demonstrate the procedure for locating the saddle point. In the following 

23 In this way one can confirm that indeed w\ — for the very symmetrical configuration described above. 




(95) 



#12 #34 
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section we will discuss the issue of the multiple saddle points and connect with the discussion 
of section 3.3 about the different configurations for the string embeddings in AdS%- 

6.3 Saddle points and multiple string configurations 

In section 2 and appendix C, we discussed the different AdS 2 string embedding geometries 
and its connection with the different boundary conditions (11) that one can impose on 7. 
At the level of the functional equations, we have seen that the different boundary conditions 
manifest in the different £ — > 0, 00 asymptotics of the coordinates. More precisely they will 
affect the constant in the expression (40). One may ask which of the configurations 
in figure 2 we should find given a cross ratio and a set of conformal dimensions. We have 
already introduced this question but let us make it more precise now that we have all the 
tools in hand. 

Consider the example of the WKB triangulation we have been studying. There exist a 
priori two choices for the orientations of the spikes, as discussed in section 3.3 and appendix 
C. Consider the spike configuration of figure 4A, which is the one used for the saddle-point 
analysis of section 6.2. Recall that for this case the cross ratio u = X2i{9 = 0) ~ 0.9 is 
positive. This means that the point 24 is located between X3 and x\. Furthermore, we have 
found a saddle point located in — 1 < w% < +1. Looking at figure 1 we see that this sit- 
uation corresponds to a case in which the insertions do not cross since x 3 < x 4 < x 2 and 
W3 < wl < wi. M Since the insertions do not cross, we expect the string embedding to be 
that of figure 2A. 

Now, one might expect that the configuration in figure 2B can be described simply by con- 
sidering the spike configuration of figure 4B but keeping the same saddle point W3 < w* A < w\ 
and the same value of the cross ratio as in the previous spike configuration. However, this 
is not possible since in this new spike configuration the cycle for X24 connects two spikes of 
the same type and thus from (40) we see that it will acquire an overall factor of (—1) so that 
u < 0. 25 

The point is that if we fix a cross ratio u and then consider a specific saddle-point w\ 
then the orientation of the spikes is fixed and thus the configuration of the string embedding 
is also fixed. This is in perfect agreement with the mapping between figure 1 and 2 and it is 
non-trivial that the integral equations encode this mapping. 

Given the above discussion, we are confronted with a very interesting possibility. Generi- 
cally we do not expect the saddle point w\ to be unique and it's likely that there are actually 
several saddle points w$ on the 4-punctured sphere. As per the above discussion, for fixed 
u, A a and a given w^) the corresponding string embedding is fixed. 26 In particular, if there 
is a saddle point in each of the three intervals of the real axis we should examine each of 

24 Be aware that the ordering of the x a in this discussion is different from that used in figure 1. 

25 Of course the corrections to X24 from iterating the integral equations will differ for the two different 
spike configurations but they should not change the overall sign of X2I ■ We are taking this as a physically 
motivated assumption in this discussion. We have checked this assumption in a few examples and found that 
it holds. 

26 Note that for fixed u, A a the triangulation will depend on which w^) one is considering. This is not a 
problem as one can apply the method of section 4 to each of these triangulations individually. 
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these in turn. 27 For a fixed u, two of these will be double-folded and one will be single-folded. 
One should find all of these saddle points and determine which is the dominant one, which 
is equivalent to ask which string embedding configuration is dominant. One may even find 
a dependence of the dominant string configuration on the dimensions of the operators thus 
giving rise to phase transitions between configurations. The issue of finding the different 
saddle points and their dependence on the parameters of the theory certainly deserves a 
deeper study. 



6.4 Extremal Limit 

In this section, we study the correlation function 

(tr Z A ( Xl ) tr Z A2 (x 2 ) tr Z A (x 3 ) tr Z Aa (x a )) (96) 
in the extremal limit when 

A 2 = 2A + A 4 . (97) 

Such correlator is protected from quantum corrections as conjectured in [16] and later proved 
in [17]. Thus, we expect to obtain the tree level gauge theory result which in the planar limit 
is simply given by Wick contractions 



™2A ™2A ~,2A 4 ' 
x 12 x 23 d/ 24 

The AdS part of our formula is universal in the sense that it only depends on the dimensions 
of the operators. On the other hand, the sphere part of the correlation function involves the 
precise details of the operators inserted. Compared to the previous sphere calculation (88), 
computing (96) just amounts to take the complex conjugate of the wave function located at 
xi, due to the replacement of Z — > Z. 

Let us start by studying the case when the cross ratio is u = 1, where we know the saddle 
point is wl = 0. From this we will be able to see the general mechanism that gives the 
expected simplification of our result. The first important observation is that in this limit 
the zeros of T(w) collide on the real axis as depicted in figure 18. Let us start by analyzing 
what this implies at the level of the x-syst em - As the integrals Uu and W34 vanish, the 
X's associated to these cycles, namely X34 an d Xi4i tend to —1. This observation has the 
remarkable consequence that the right hand side of all equations in the x-system becomes 
trivially equal to 1 as one can easily verify 28 . As a result, all x-functions are exactly given by 

27 Here we are not considering the possibility of complex-conjugate pairs of saddle points located off of the 
real axis (see footnote 2). 

28 This trivialization of the x-system is general and follows just from the fact that the two cycles W14 and 
W34 vanish which implies that the %-functions X34 an d Xi4 become -1. In the specific case of U = w 4 = 
and Ai = A3, which turns out to correspond to cross ratio 1, the x-system is already trivial because of 
the symmetry of the stress energy tensor in this particular point of the parameter space, see footnote 16). 
Nevertheless, we emphasize that the trivialization of the x-system in general does not rely on this specific 
symmetry of the stress energy tensor. 
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Figure 18: In the extremal limit, the main feature is that the zeros collide on the real axis. The 
black lines represent the WKB cells whereas the gray lines represent the WKB triangulation. At 
the exact extremal configuration, there are no WKB lines connecting 1 to nor —1 to 0. We 
interpret this as a manifestation of the field theory fact that at tree level all operators are Wick 
contracted only to the fourth operator. 



leading term of the WKB expansion (40) 29 . For convenience, let us introduce an infinitesimal 
S denned by the condition 5 = 2A + A 4 — A 2 . At the end of the day, we will take S — > 0. In 
this limit, the solutions of the x-system are then given by 

X23 — X12 — -e 2 , X34 — Xi4 — -e 2 , X24 — X 24 — e l yy J 

We may now plug this solution in the expression (51) and extract the cycles using as described 
in section 4.10. We find that all A's vanish in the limit 5 — > except for v4 14 and A 34 , which 
tend to —1 as S goes to zero. This implies that all rjE ab vanish except for 7714 and 7/3,4, which 
diverge since the integrand of these cycles becomes singular in this limit. However, one must 
go back to the area formula (66) and realize that such cycles are multiplied by a vanishing 
quantity. Indeed, (66) simplifies to 

^ 14 + |. (100) 
In the limit S — > 0, the first term of this expression is explicitly given by 

d6 cosh 9 log (l - e -|^ coshe J + 0{5) = — - + 0{5). (101) 

Hence, it turns out that the finite AdS contribution vanishes in the extremal limit. We 
believe this is the general mechanism for any value of the cross ratio. 

The computation of the sphere contribution follows the same steps as before, with a slight 
change on one vertex operator (recall that to get the extremal case, we replaced the operator 
located at x\ in (88) by tr Z A ). The new solution for the equations of motion is 

(p(w,w) = i(A\og\w + 1| + A 4 logH + Alog|u7 - 1|) . (102) 



29 Indeed, when the right hand side of the %-system is 1, the kernel term in equation (54) vanishes and we 
are left with leading WKB contribution. 
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Now when we compute the contribution of the sphere action and wavefunctions on this 
solution, we find that it exactly cancels the term y/TT for A's satisfying (97). Consequently, 
the sphere part of the correlation function also vanishes in the extremal limit. 

The divergent piece in the extremal becomes simply 

e-^° f[(\s a A s b \ )-^ -+ 5 log(l-e^ c °^)+0(5) = log ^+0(5) 

(103) 

which goes to zero as 5 — > 0. We are left with the spatial dependent part which, using that 
the cross-ratio is 1, can be written as 



(^2) 2A {^f)^ {^f) 2 ^' 

This is nothing but the tree level result (98) of the gauge theory. 



7 Discussion and future directions 

In this paper, we have computed the AdS part of the four point function for heavy scalar 
operators in M = 4 SYM in the classical limit. For the particular case of BPS operators on 
a line with a single scalar field, the sphere part is known and thus we can construct the full 
strong coupling four point function. 

The main ingredient of our method is the integrability of the string equations of motion 
in AdS2- Specifically, we use the method of Pohlmeyer reduction to map the problem to that 
of solving a certain modified Sinh-Gordon equation which is known to be integrable. We 
construct the linear problem associated with this equation, which has the form of an SU (2) 
Hitchin system. This approach was used in the solution of the Null Polygonal Wilson-Loop 
problem at strong coupling [10] as well as in the study of three-point functions of heavy 
operators at strong coupling [8, 9, 11]. 

Let us mention that while our approach was inspired by these previous works, to solve 
the N > 3 point function problem required significant generalization of [10, 8, 9] as well as 
nontrivial new ingredients. For the case of the Null-Polygonal Wilson loop the world-sheet 
has the topology of a disk, whereas in our problem it is that of an iV-punctured sphere and 
this changes the boundary conditions that one imposes. This issue was addressed in [8, 9] for 
the case of the 3-punctured sphere, however in those works the total monodromy condition 
was enough to derive the functional equations that determine the necessary objects. These 
functional equations are linear and can be easily inverted using standard techniques. For the 
case of 4 or more punctures the situation becomes significantly more complex. First of all, 
the total monodromy condition no longer provides enough information to fix the necessary 
objects. We have made heavy use of the formalism developed in [12] to derive the functional 
equations. Second, the inversion of these functional equations is more subtle due to their 
complexity. The result turns out to be some integral equations resembling the usual TBA 
equations. 
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Perhaps the most exciting aspect of this work is the multitude of interesting applications 
and extensions of the results. Let us consider each of these in turn. 

• Multiple configurations and phase transitions. An important physical outcome of this 
paper is the emergence of multiple string configurations in AdS2- Each of these config- 
urations is associated to the existence of several saddle points. A natural question is to 
figure out whether the dominant saddle point depends on the parameters of the theory. 
If so it would be interesting to study the phase diagram and the possible transitions. 
We have already made some preliminary progress in this direction and we hope to make 
a more extensive study in a future publication. 

• OPE. A natural question to ask given any 4-point function in a conformal field theory is 
what can be learned from its OPE decomposition. In particular, important information 
about the spectrum and structure constants of the theory can be extracted. 

• GKP string. An interesting aspect of [8, 9] is the similarity between the mathematical 
formalism employed despite the differences in the physical problem: [8] describes strings 
without spin in AdS2 whereas [9] describes spinning strings in AdS^. In the formalism 
of [9, 11], one expresses the iV-point function of GKP string in terms of a universal AdS 
contribution and a contribution from vertex operators, both of which can be computed 
for the case of the three-point function. It is possible that one could use the formalism 
developed in this work to calculate the AdS contribution to the iV-point function of 
GKP strings. 

• N -point functions. The formalism developed here does not depend in any special 
way on having only 4 punctures and in principle one could use the same methods to 
study the iV-punctured sphere for any N. It would be interesting to understand how 
the functional equations generalize to higher N . Furthermore, since for the iV-point 
function there will be N — 3 unfixed insertion points, the moduli space of possible 
configurations should be quite interesting. 

• TBA equations. We should note that the techniques developed in [12], in principle, 
allow one to write the functional equations derived in this paper in the usual form of 
a Y-system. Typically this Y-system will involve an infinite number of Y-functions. 
This form of the equations could be useful for various applications including analytic 
continuation of parameters and generalization to iV-point functions. 

• Generalizing out of the line and WL/CF duality A natural step would be to generalize 
this work for operators not inserted on a line. In this case the string is embedded in a 
higher dimensional AdS space, which involves a more complicated Pohlmeyer reduction 
scheme. It would be interesting to study the question of whether the multiple string 
configurations/ saddle points we have found is special to AdS% case. Another promising 
application of such generalization would be the possibility of studying the OPE for the 
Null Polygonal Wilson Loop [22]- [24]. One could also investigate the duality between 
Null Polygonal Wilson Loops and Correlation functions of null separated local operators 
at strong coupling [25]- [30]. 
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Many of these points present interesting opportunities to try to learn about finite coupling 
and weak/strong coupling interpolation and this is probably the most stimulating reason for 
pursuing them. 
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A The linear problem 

A.l Summary of the linear problem 

The linear problem associated with (10) is given by 



(d + J w ) i> = 0, (8 + J m ) ip = 



(105) 



where the connection has the form 




w 



&w + A. 



w 



(106) 



w 




(107) 




(108) 
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w 




(109) 




(110) 
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For compactness we have introduced the combination 7 = 1/2(7 + log V TT ) . The function 
7 is defined as the solution of the following problem 

dd'y = \j TT sinh 7 

7 -> i^logTf {w->Z a ) (111) 
7^0 (tO->!O a ) 

where z a and w a are the zeros and the poles of T, respectively. 

For the near-puncture analysis as well as the WKB analysis it is useful to make the field 
redefinition ip — >■ ip — Gip where 

A _ 1 ( +e -7/2 T l/4y-l/4 1 \ , . 

L '" <1 _ < ,-r/2 T l/4*-l/4 1 I l 11 ^ 



2 y _ e -7/2yl/4 T -l/4 2 

This is usually refered to as 'diagonal gauge' in the literature. In diagonal gauge we have 

= ^(-"o 1 ?) (113) 
1 

2 v x V — sinh 7 cosh 7 . 



<§>_ - -\/t ( - cosh 'y sinh T ^ fi 14^ 

— o V V -sinh 7 cosh7 ) \ 1L ^J 



Aw ~ 9w ( 47 -|t (TT) (115) 

A ™ - <M i 7 +ilog(Tf)J ^ iib J 



We are now ready to consider the behavior of the solutions near the points w a and z a . 



A. 2 Solutions near w a 

Let us first consider the solutions of the linear problem in the neighborhood of one of the 
punctures. From (111) and the explicit expressions for <3> and A for w — > w a we have 

^-►^(VA). Sa-^v^UA) (ii7) 

Aw -> ^ (~ logTf ) ( + 1 £ ) , ,5* ^ Q log Tt ) ( V £ ) (118) 
Then the solution in the vicinity of puncture P a is given by: 

^ H = (t/t ) 1/8 e ±irr^| ±) „ (w _ ^^lA.r 1 -! (u7 _ ^ a fl^+l \ ±) ( n9 ) 

where |±) are the eigenvectors of the Pauli matrix a 3 . Note the characteristic monodromy 
of the solutions about w a . 
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A. 3 Solutions near z 



Now we will consider the behavior of the solutions near the zeros z a of T. Notice from (107) 
- (110) and (111) that the connection is regular or singular at z a depending on the direction 
of the spike in 7 at z a . More specifically, the connection is regular if 7 ~ — log \ w — z a \ and 
thus the solution will be regular in the vicinity of a <i-spike. However, the connection has a 
singularity if 7 ~ + log \w — z a \ and at the w-spikes one can check that in gauge (107) - (110) 
there are two linearly independent solutions behaving as 

* ~ **» = I ( W - Za )^\*-- Za )-^ ) (120) 

where we have written the two solutions in matrix form as in (71). Notice that \I> has square- 
root type singularity at z a since it has a monodromy of \I/ — > (—1) ^ about z a . The solutions 
associated with the punctures {sp} and {Sp} inherit this square-root singularity as one can 
see by expanding them in the basis (120) near z a . 

In our analysis it is crucial to account for the additional monodromies originating from u- 
spikes. Let us explain our conventions for doing this. If there is a w-spike at z a , one can always 
make the gauge-transformation \1/ — > y!*" 1 ^ that removes the square root singularity (\l/ 2a 
is given in (120)). Of course this gauge transformation contains the same multivaluedness 
and one must still account for it at the end of the day. In the main text we use the point of 
view that this gauge transformation has been performed for each u-spike. The connection in 
this gauge will only have singularities at the punctures and the solutions in this gauge will 
only have non-trivial monodromies around the punctures. In this way we can define small 
solutions that are single valued throughout some Qe, as is the prescription of [12]. We must 
then be sure to account for the multivaluedness of these gauge transformations whenever 
we have a holonomy that encloses an odd number of w-spikes. Such holonomies arise in the 
WKB expansion of the coordinates and we will return to this issue below. 



B WKB analysis 

B.l Statement of the WKB approximation 

As we have discussed above, it is essential to have control over the £ — > 0, 00 asymptotics 
of the inner products. It is clear from (105 — 106) that these are both singular limits, and 
the basic idea of extracting this singularity is as follows. As discussed above, we have good 
control over the solutions in the neighborhood of the punctures. Thus what we would like 
to study is the transport 



Pexp 



$ + A + 



i> K) 



(121) 



where C {w' a — > w) is a curve starting at w' a , a point in the neighborhood of w a , and termi- 
nating at a generic point w. Note that at any point on the punctured sphere C the Higgs 
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field $ has the two eigenvalues =Fw/2 = =pvT'/2 dw (which are single valued on the double 
cover E), and thus we can choose a gauge along C where $ is diagonal and given by 

®=-(~n °) (122) 

Now consider the £ — > limit. First consider an infinitesimal segment of C in the neighbor- 
hood of P a . In the neighborhood of P a the connection (in diagonal gauge) becomes diagonal 
(see (117)-(118)) and thus one can break apart the path-ordered exponential. In particular, 
one can isolate the singular part e - ^ */^|±) which will have one component growing expo- 
nentially and one component decaying. Let us choose the branch of $ such that the |+) 
component is the one that is growing as we transport along C away from P a (although for 
the moment we are still working in a neighborhood of P a ). This will correspond to the small 
solution at P a since it is exponentially decaying as it is transported toward from P a . The 
WKB approximation is the statement that the exponentially growing part of the solution as 
£ — > will continue to be given by e~f®^\+} as we transport away from the neighborhood 
of P a (now leaving the neighborhood of P a ) as long as we follow a curve such that at every 
point we have 

Re > (123) 

This condition is satisfied most strongly along a curve such that 

Im(w/O=0 (124) 

Condition (123) is called the WKB condition and curves satisfying (124) are called WKB 
curves [12]. Along a WKB curve defined for Arg(£) = 9 the WKB condition is satisfied 
for Arg(£) £ (9 — tt/2,6 + 7r/2) and the WKB approximation is guaranteed to hold in this 
range. For example, suppose there is a WKB line connecting P a to Pj, for 9 £ (6L,# + ) 
but not outside that range. Then the WKB approximation will reliably give the £ — > 0, oo 
asymptotic for 9 £ (0_ — tt/2,9 + + n/2). These statements are proven in [12] and we refer 
the reader there for a more detailed discussion. 



B.2 Subleading WKB 

We will now consider the £ — > limit of the inner products (or Wronskians) (s& A s a ) (£). We 
consider the case when P a and Pj, are connected by a WKB line which will be an edge E a b 
in the WKB triangulation. From the analysis of A. 2 we know s a and s& in the neighborhood 
of P a and Pb respectively. In order to evaluate the Wronskian we need to know the solutions 
at a common point. The approach here is to use the connection to transport the solution s a 
along E a i) to a point w' b in the neighborhood of P& and then to evaluate the Wronskian at 
w' b . That is, we want to study the £ — > behavior of 



(s 6 |Pexp 



dt -H + V 
o 4 



(125) 



where we defined 



H = w$ l 



V = wA w + wAyr, + £w$, 



(126) 
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The contour of integration in (125) is the edge E a b and the components of (126) are defined 
in appendix A. The basic idea of the computation is to expand in a perturbative series where 
^~ X H§ acts as the free Hamiltonian. Such a procedure will be valid so long as the free part of 
the Hamiltonian is sufficiently larger than V for all points along the curve, which will be true 
along the edges of the WKB triangulation. Then we can expand (125) in the Born series 



fo H o/t\ 



+} 



Jo 



(127) 



+ ! dt 2 ! dhi+le-X 2 Ho ^V{t 2 )e'^ Ho ^V{t x )e-^ Ho/ ^\+) 



Let us explain a subtle point regarding the 'external states' in the above expression. We start 
with the small solution at P a which we take to be We then transport it to Pb and then 
extract the coefficient of the exponentially growing part - that is, we take the inner product 
with the small part of this transported solution. Since ~ 1+) grows as we transport it 
along a WKB curve (i.e. it decays as one follows the curve into P a and thus grows as we 
transport it away from P a ) and Ho is diagonal, we infer that the small part of the solution 
at Pb is the solution proportional to |— ). Thus we take the out-state to be (— . Finally, 
since the inner product is the antisymmetric the (— | gets flipped to a (+|. 

Using the fact that |±) are eigenstates of the free Hamiltonian we can easily evaluate 
the order O(V ) and OiV 1 ) terms in (127). For the 0(V 2 ) term, we insert the identity 
+ I between the two insertions of V. We find 



(-l)^e + i/o^| 



1 



dt 1 (+\V(t 1 )\+) + 



o 



dti(+\V(ti)\+) 



+ 



[128) 



i ft 2 
dt 2 I 
o Jo 



dt 1 e-%^(+\V(h)\-)(-\V(t 1 )\+) 



Now concentrate on the second term on the O (V 2 ) contribution. As £ — > the factor 

exp ^— J ( * 2 uj/^j will suppress the integrand except for the small range t 2 = ti + 0(£) and 

thus the result of the first integration will already be 0(£). So to order £ we can take u to 
be constant and V(t\) — > V(t 2 ). We then find for the second term in the 0(V 2 ) contribution 

fftfMM (129) 
Jo wy/T 

Putting everything together, we see that the result re-exponentiates and we find 



'~ l )^b eX P 



+ -J dt-VT- I dt{+\V(t)\+}+£ 



1 dt \(Mv(t)\-W 



wyT 



(130) 



Grouping each term based on its order in £ (including the prefactors ip b whose explicit 
expression are given in (119)) we find 



{sb A s a ) (f ) ~ exp 



^ 1 ZU ab + a ab + ^TUab + iVab 



(131) 
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where 



<y a b 



Vab 



lim lim 

w' a ->w a w' b -+w b 



I Vfdw + \og(w a -w' a ) + ^Y hg(w b - w' b/ 
I ( \ dw ( 7 ~ log dw + \ d *>(l + log Vtt) dw) 



1 r= 1 2 

-VT(cosli7 — 1) g?u> H —(dj) dw 

2 4VT 



(132) 
(133) 
(134) 



This completes the derivation of formula (26) used in the main text. The integral w a b is 
defined as in (58). Note that the logarithmic terms in w a b in (131) are due to the prefactor 
■0k~V>a ■ These terms precisely cancel the singularities at the endpoints of the integral J^u so 
that w a b is finite as we continue the limits of integration all the way up the to punctures at 
w a and Wb [8]. 



B.3 WKB expansion of the coordinates 

In the previous section we derived the £ — > WKB expansion of (s a A s&) up to order 
To compute the WKB expansion of the coordinate xe we simply combine the expansions 
for each edge of the quadrilateral Qe, taking care to account for the directions of the WKB 
lines as discussed in section 4.6. When this is done each of the integrals (132) - (134) become 
closed integrals along the cycle je- The asymptotics of the x-f unc tio ns are needed for the 
inversion of the x-system described in section 4.9. For that purpose only the non-vanishing 
contributions are needed in the £ — > 0, oo limits. 

There is one very important subtlety that must be addressed here, which is that of the 
monodromy around w-spikes discussed in appendix A. 3. We take the point of view that 
we have made the (multi- valued) gauge-transformation (120) that removes the monodromy 
about each w-spike. The small solutions in this gauge are single valued throughout Qe, but we 
must account for the monodromy of the gauge transformation about Qe- This monodromy 
is simply (— 1) UE where ue is the number of u-spikes in Qe- 

Combining the above discussion with (133), the constant term in the WKB expansion of 
Xe is given by 

Cf = log(-l) UB -j J (dj+ *dlog Vrf) = log(-l) MB ± in (135) 

To arrive at the last equality (135) we have used the fact that 7 is single- valued on the 
4-punctured sphere so that the integral of dj on any closed contour is zero. The integral of 
*d\og VTT is simple to do explicitly and gives the ±in factor. 30 

The discussion of the £ — > 00 limit follows along the same lines as the £ — >■ limit. 

30 The ± depends on the orientation of je but both signs have the same overall effect so that the ± is 
irrelevant. 
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The singular term is given by e^s^ 2 . 



The constant term is the same. Thus the full 



non-vanishing WKB asymptotic is given by 

X E~(-l) UE exp 



If (rw^) 

Z J IE 



(136) 



where we recall that ue is the number of w-spikes enclosed in je- This is the expression (40) 
used in the main text. 



C Fold lines and Properties of 7 

In this appendix we discuss some properties of the function 7 and how they are related to 
geometric features of the string embedding. In appendix C.l we show that the world-sheet 
contours where 7 = map to the fold-lines of the target space solution; in appendix C.2 we 
discuss how the geometry of the string embedding near the boundary is deduced from the 
structure of these 7 = contours near the points w a ; finally, in section C.3 we show how 
the global structure of the 7 = contours is deduced from the choice of spikes in 7. The 
point of this appendix is to give the background details that were omitted in the discussion 
of section 3.3. 



C.l Fold lines 

In this section we show that the contours on the worldsheet where 7 = map to the fold 
lines of the string embedding. This was pointed out in [8]. Recall the relation between 7 
and the world-sheet metric 

/—= , dxdx + dzdz , „. 

V TT cosh 7 = 137 

z A 

Furthermore, we have 

(«.)■+ (a,)' m= m*+w {m 

z 2 z 2 
Now, suppose that C is a curve on the worldsheet that maps to a fold-line of the string and 
consider a point V in that curve. We can choose local coordinates at V so that the derivative 
takes the form d w — > e 1 ^ (d t + id n ) where the direction d n is chosen such that d n z = 0. 
The prefactor is the Jacobian of the coordinate transformation (just a translation and 
rotation). The defining property of the fold-line is then that the x-coordinate reaches an 
local extrema and thus we also have d n x = as we cross the fold. Therefore along the 
fold-line we have (with dtx = x) 

1 /-y»2 I 2 /y>2 I ^2 ,-yi 2 I ^,2 



VTT cosh 7 -)■ T(w) -)• e^—^, f(w) -> e -^-^JL_ (139) 

z l z l z 1 

Using the last two equations to solve for VTT we see that they are consistent with the first 
equation only if 7 = 0. Therefore, the worldsheet contours where 7 = map to the fold-lines 
of the string-embedding. For this reason, we frequently refer to the contours where 7 = as 
fold-lines. 
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A MB 

Figure 19: Single-folded and double-folded string in panels A and B respectively. 

C.2 Structure of 7 near w a 

To gain some intuition about the structure of the contours where 7 = it is useful to study 
the behavior of 7 near the points w a . Recall that 7 is defined as the solution to the PDE: 

957=v / TTsinh7 (140) 

subject to the boundary conditions 

7 -> i^logTf (w^z a ) (141) 

7^0 (w^P a ) (142) 

The boundary condition (142) simply imposes that 7 is non-singular at the singularities of 
T and this condition is automatically imposed if we demand the solution be regular away 
from the zeros of T. 

Since we know that 7 must vanish at singularities of T, it's natural to study the function 
in the neighborhood of these points. Let us consider some P a and use polar coordinates 
(r, 0) in which the origin is at w a . Since 7 is vanishing, we can linearize the RHS of (140). 
Further, we can take VTT ~ |A| 2 /(4r 2 ). The PDE becomes linear and separable and using 
standard techniques one finds the series solution 

00 

7 ~ 9o r^ A + 9m sin (m0 + S m ) r^ AW (143) 

m=l 

Now consider a small circle centered at r = 0. As r — > the series (143) is dominated by the 
lowest mode in the expansion. Thus along an infinitesimal circle centered at r = the series 
(143) will vanish 2m* times, where g m * is the smallest non-zero coefficient g m , m — 0, 1, 2, ... 
in the series. Thus, if go is the smallest non-zero coefficient then the series will vanish only 
at the point w a which will be a local extrema. If m* = 1 then the series will vanish along a 
single curve passing through P a ; if m* = 2 then 7 will vanish along two curves that intersect 
at P a , and so on. 

The fact that the contours where 7 = map to fold-lines of the target-space solution 
gives a clear geometric meaning to each possible behavior m* = 1,2,... near an insertion 
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point. For m* = 1 we will cross two fold-lines as the world-sheet coordinate traverses a small 
loop around the point w a . This means that near the insertion point the string is single-folded 
as shown in figure 19 A. For m* = 2 we will cross two fold-lines as the world-sheet coordinate 
traverses a small loop around the point w a . This means that near the insertion point the 
string is double-folded as shown in figure 195, for example. In general for n > the case 
m* = n should correspond to an n-folded string. The only subtle case seems to be m* = 0. 
Apparently if m* = 0, as we traverse a closed loop around w a the contour swept out in the 
target space does not close since there is no point at which the coordinates (x, z) can 'turn 
around'. In this paper we are only interested in solutions that are closed (i.e. the embedding 
coordinates have trivial monodromies around operator the insertion points x a ) and thus we 
will only study cases for which m* > at all w a . This is further discussed in appendix C.3. 

It is important to keep in mind that (as we mentioned above) the behavior of 7 at P a is 
not our choice, and is determined by regularity and the conditions (141). In other words, for 
fixed T the only remaining conditions one can specify are the choice of signs in (141). For 
each choice of signs there will be a unique m* for each P a . In the next section we demonstrate 
how this works using the T of the 4-point function discussed in the main text. 



C.3 Structure of contours where 7 = 

In this section we describe why the spike configurations of figure 4 are the only two physically 
relevant configurations. Furthermore, we deduce the structure of the contours where 7 = 
for each of these spike configurations. 

Consider T fixed to be that of the 4-point function discussed in the main text (see equation 
(60) and figure 15). There are 4 zeros and therefore 2 4 ways to choose the signs in (141). 
Because of the symmetry of (140) under 7 — > —7, without loss of generality we can fix one of 
the spikes to be up which leaves 2 3 choices. Now, because the string is embedded in AdS2 we 
know that it must be folded. Moreover we know that the operator insertions x a will sit along 
the fold-lines of the target-space solution. In the world-sheet coordinates this translates to 
the statement that we should require m* > at each w a . That is, there should be at least 
one contour where 7 = running through each insertion point w a . For the 4-point function 
T (see equation (60) and figure 15) the only obvious way to accomplish this in general is to 
choose the spikes such that 7 — > —7 under reflection about the real axis. This leaves only the 
spike configurations of figure 20A,B, which are those of figure 4 used in the main text. We 
will now discuss the global structure of the 7 = contours for these two choices of spikes. 

In figure 20 we show the fold-structure for three different spike configurations. The 
black lines schematically represent the contours where 7 = and one can read off the m* 
associated with each puncture. The structure of these contours is determined purely by the 
choice of the directions of the spikes of 7. We refer to these contours as 'fold lines' since they 
map onto the fold-lines of the target-space embedding (see appendix C.l). We guess the 
structure of the fold lines for each choice of the spikes as follows: u spikes must be separated 
from d spikes by at least one fold line; we use the minimum number of fold lines needed to 
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m* = 1 m* = 1 m* = 1 m* = 1 

X 



d 
d 




Figure 20: Three different spike configurations and the corresponding structure of the 7 = 
contours. The black lines schematically represent the contours where 7 = and one can read off 
the m* associated with each puncture. Panels A and B show the physically relevant configurations 
studied in the main text. Panel C shows a third spike configuration which is not physical due to 
the presence of m* = behavior at two of the insertion points. The gray contours in panel C 
indicate contours that cannot cannot correspond to fold-lines due to the restriction that 7 = 
contours must encircle at least one zero of T (see footnote 31). In this figure we are not indicating 
the location saddle-point w\ because it is not relevant for the present discussion (so long as it is 
located somewhere on the real axis). 



accomplish this for all spikes. Note that fold lines must encircle at least one zero of T. 31 This 
restriction is useful because, for example, it allows one to rule-out the possibility of fold-lines 
corresponding to the gray contours in figure 20C. This is important because if it was possible 
for the gray contours to be fold-lines then it might be possible to have a solutions with all 
m* > for configuration C. Configurations A and B are the physical configurations that 

31 Consider a closed contour along which 7 = and suppose (for a contradiction) that it does not enclose 
any zeros of T. Let T> be the region enclosed by the contour. This contour must separate positive values 
from negative values (i.e. it cannot sit at the bottom of a 'valley' since this locally violates the equation 
(140)). Suppose for simplicity that 7 < in T>. Since 7 is regular away from the zeros of T, there must be 
at least one local minimum inside T>, and therefore at least one point where (d%. + dy) 7 > 0. Thus at such 
a point the LHS of (140) is positive or zero, but the RHS is strictly less than zero by assumption, which is 
the desired contradiction. 
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we study in this paper and we have checked the fold structures of figure 20A,B numerically. 
Configuration C is an example of a spike-configuration that does not correspond to a target- 
space solution with the desired properties; the corresponding fold structure is only our best 
guess but we have not checked it numerically 

To summarize this appendix, in appendix C.l we showed that the world-sheet contours 
where 7 = map to the fold-lines of the target space solution; in appendix C.2 we discussed 
how the geometry of the string embedding near the boundary is deduced from the structure 
of these 7 = contours near the points w a ; finally in section C.3 we discussed how the global 
structure of the 7 = contours is deduced from the choice of spikes in 7. From all of this 
one can deduce some qualitative global features of the string embedding, which is discussed 
in detail in section 3.3. 



D Details of the 4-point function computation 

D.l Explicit expression for stress-energy tensor coefficients 

For completeness we present the coefficients c a of the stress-energy tensor in formula (60), 

A 2 

00 



4 



c = - [4:Uw4 + 2w4 (1 + w 4 ) A 2 + (-1 + wa) (2-u; 4 A 2 + (1 + w 4 ) (A 2 - A 4 ))] 
ci = ^[-2f/+(-l + W4 ) 2 A 2 -(l + W4 ) 2 A 2 ] 

c 2 = X - [-4£/w 4 + 2(l + u; 4 )A 2 + (-l + w 4 ) (-2A 2 + (1 + w 4 ) (-A 2 + A 2 ))] (144) 
D.2 Explicit expressions for ^-functions and Apq 

For reference, we include here the explicit expressions for the x-functions for the triangulation 
of figure 15. They are given by 

(si A Mf 1 g 4 )(g 2 A S4) 
(M 1 s 4 A S2){s 4: A si) 

X23 = (-l) (g2AM3g4)(g3AS4) (146) 

, ^ (s 4 As 2 ) (s 3 AM^s 2 ) 
(s 2 A s 3 ) (M 3 s 2 A s 4 ) 

, (s 4 AM lS2 )( Sl As 2 ) 
Xl4 " { - 1) (M 1 s 2 As 1 )(s 2 As 4 ) (148) 
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(S 2 As 3 )(s A A St) 
X2A = ( — IJ7 77 7 (149) 

(M^s 2 AM ASl )( Si As 3 ) 

One can check that these coordinates satisfy the rule (37) at each puncture. The v- s yst em 
obeyed by these coordinates is given by 

v v++ fv - Y++V 1 - (I±M1I±M nsn 



;i + a 34 ) (1 + a 12 ) 

1 + .1, 

^24; 



X12X++ = (Xl4X^)" 1 = X34Xi + =(x 2 3X2 f 3 + ) _1 = |^ 1 ( (152) 



iven 


by 




Xl2 ( 


;i + xu) 


(1 + X24 (1 + X23 (1 + X24))) 


X23 ( 


(i- M f) (i - 
;i + X34) (1 + X24 (1 + X12 (1 + X24))) 


X34 ( 


'1 + X23) 


(i-mI) (i - mI) 

(1 + X2 4 (1 + V 14 (l + X24))) 


Xl4 ( 


;i + xi 2 ) 


(i - mI) (i - Ml) 

(1 + X 24 (1 + X34 (1 + X24))) 



. A12 V x ' XU) \ x ' A24 V x ' A.23 V T X24JJJ 

(1- Ml) (1-^2) 
, _ X23 (1 + X34) (1 + X24 (1 + Xl2 (1 + X24))) 

A23 ~ (l-M^(l-Ml) (154) 

. X34 (1 + X23) (1 + X24 (1 + Xl4 (1 + X24))) 

^34 = 77 2T7^ 2\ ( 155 ) 



1 All \" 1 Al^y V 1 A,<a« 1 AOI V- 1 1 /\A1J J J r«\ 

A 14 — 77 2\ [LOO 

A X24 (1 + Xl2 (1 + X24 (! + X23))) (1 + V 43 (1 + A42 (! + Xu))) n K ^ 

Mi - (HIM (157) 

. X24 (1 + X23 (1 + X24 (1 + X12))) (1 + Xu (1 + X42 (1 + X43))) ( . K o\ 

(i-mD(i-mI) ( ] 

Using the explicit expressions for the coordinates (145)-(150), schouten identity and the shift 
relation (43) one can directly verify the functional equations (151)-(152). 



D.3 Finite part of AdS 

In this section, we present some intermediate steps in the derivation of our formula (66) for 
the finite part of the AdS contribution. We want to compute 




(159) 



according to the steps outline in section 5.1.3. The complete basis of five a-cycles and five b- 
cycles that we chose is depicted in figure 21. From this figure we also read-off the intersection 
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Figure 21: The cycles for Riemann bilinear identity. The dashed line represents a contour in a 
different Riemann sheet. The wavy lines represent a choice of branch cuts. From this picture 
we also read the intersection matrix I a b of the cycles. For each pair of cycles, say 7 a and 75, 
intersecting at a point with tangent vectors d a and <9 respectively, we assign I ab = +1 (—1) if 
det [{d a ,d b }} >0(<0). 



matrix I ab = (5 a +i,& — #a-i,b) using the conventions described in the caption. The only other 
ingredient we need is 

,7 = 0, z = 2,...,5 (160) 

which follows from the regularity of 77 at the poles of T. Plugging into (159) and computing 
we find 



.4 
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fin 



■ + i {u ai Vz 3 ,z 2 + W« 2 ^l, 22 + U a3 T) ZyZ2 + U] ai T] 1)Z2 + 0; a5 7700,. 



5^w 6j r] Z3tZ4 (161) 



,i=l 



where we are using the notation r/ a b = rj and ujq = J c u) and the contours are defined in 
figures 21 and 22. 

Each of these f] a b can be written as a linear combination of the rjE ab = j E rj where the 
integral is taken along the WKB-line from P a to P b and the direction of the contour is the 
same as that of the WKB line. The idea is to combine the r/E ab to form the contour that 
we want. Let us exemplify with rji jZ2 . From the WKB configuration, see figure 22, we see 
that the large 6 expansion of the ratio ^ lA ^^^ AS4 - ) involves a cycle that can be continuously 
deformed into twice the line integral connecting the puncture at w — 1 and the zero at 
w = Z2- Therefore we have 



~2 



de 



71 



'log 



(1 + A 2A ) 



{Ve 12 + Ve 14 ~Ve 2 



(162) 



In the same way we obtain 
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Figure 22: To extract line integrals connecting a zero to puncture or connecting two zeros we 
combine products of elementary solutions that have WKB expansions involving integrals over the 
paths indicated by the black lines. The resulting closed contours can be continuously deformed 
into the contour that we want, indicated by the green lines. The precise way of combining these 
products is dictated by the direction of the WKB lines indicated by the gray arrows. 



V-l,z 2 

Vz,z 2 

VoO,Z2 

Vz 3 ,z 2 



2 (V24 - Vid 



(7712 - Vu - ^24) 

{Vu - V12 - 

{Vn + Vu - Vu - V23) 



(163) 
(164) 
(165) 
(166) 
(167) 



where the notation is the natural simplification of that used in (162). Plugging these ex- 
pressions into (161) and re-collecting each rjE, one finds that the coefficient of T]e is simply 
ue where uje is the w-cycle the intersects edge E, not the integral of u along edge E (which 
would be divergent). That is, it's (1/2 of) the w-cycle associated with the coordinate \e 
which are shown in figure 15. Thus we have 



A 
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fin 



— I 



E 



ue Ve 



(168) 



which is formula (66) as desired. 

Equation (168) is perhaps the simplest possible result one could write from the triangu- 
lation data. Given this simplicity, it is probably possible to derive the result in a much more 
elegant way and perhaps even for any number of punctures. We have not pursued this issue 
but feel that it merits further exploration. 



E Three-point function in GMN language 

In this section we apply the method developed in section 4 to the three point correlation 
function studied in [8]. We use the setup of [8], namely the same stress-energy tensor. We 
aim at deriving a set of functional equations to extract the cycles used there. 
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P2 



P3 




Pi 



Figure 23: The WKB triangulation for the 3-point function is composed of 3 edges forming two 
triangles on the sphere. Here we show the construction of the coordinate xi3- We are using the 
edge-splitting procedure discussed in section 4.5 (in particular, see figure 12). The gray contour 
shows the cycle associated with the coordinate xi3- 



As a starting point, we introduce the WKB triangulation for this configuration from 
which we define the coordinates, see figure 23. From this figure, we easily derive the \- 
system. Since the quadrilateral is very degenerate it follows from (46) that the right hand 
side of the x-system is equal to 1. The reason is that the same auxiliary Apg's appear both 
in numerator and denominator canceling each other. Hence, the solution of the functional 
equations is simply given by the WKB asymptotics. More explicitly, the x functions take 
the form 



where a, b and c are distinct. 32 The spikes must be in pointing in opposite direction as 
follows from the discussion of appendix C. This is the the origin of the (—1) pref actor in 
(169). The cycles of u are given in terms of the dimensions of the operators, 



Having the solutions of the functional equations, we can easily find the auxiliary quantities 
Apq using the rules of section 4.8. The determination of the 77-cycles is also straightforward. 
To compare with the result in [8] let us set Ai = A 2 = A and A 3 = Aoo. We use expression 

32 This result also follows directly from the definition of the coordinates in terms of the small solutions, 
Xac = ~ rM^^T^7~ a Sb \ f° r distinct a, b and c; all the inner-products cancel and one is left with only the 




(169) 




(170) 



Aac (M- 1 s b /\s c ){s b /\s a ) 

monodromy factors in (f69). 
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57 to compute the cycles, and we get 

/l r 

V = —e- 9 log{l + A u ) = h{2A-A 00 ) + h{2A + A OD )-2h{2A) (171) 
-l Jr * 

POO P 

/ V = —e- 9 \og(l + A^) = h(A 00 ) + h{2A + A 00 )-h(2A)-h(2A 00 ) (172) 
h Jr ^ 



h(x) 



where we define 

p jn 

/ — cosh 9 log (1 - e -^ coshe ) . (173) 
Jr 71 

This is precisely the result obtained in [8] . A last comment about the expression for the area 
in the three point function. It is easy to show using the same type of manipulation of the 
four point function case that the area can be expressed in terms of elements of the WKB 
triangulation as 

A fin = g ~ i ^ veVe (174) 
e&t 

where the sum is over the edges of the triangulation of figure 23. As in the case of the four 
point function, we define i]E ab as the 77-cycle that passes along edge E a f, from P a to and 
u>E ab as the cu-cycle that intersects edge E ab . 



References 

[1] J. M. Maldacena, "The Large N limit of superconformal field theories and supergravity," 
Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] [arXiv:hep- 
th/9711200]. 

[2] J. Escobedo, N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Functions and 
Integrability II. Weak/strong coupling match," JHEP 1109, 029 (2011) [arXiv:1104.5501 
[hep-th]]. 

[3] J. Caetano and J. Escobedo, "On four-point functions and integrability in N=4 SYM: 
from weak to strong coupling," JHEP 1109, 080 (2011) [arXiv: 1107.5580 [hep-th]]. 

[4] N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Functions and Integrability 
III. Classical Tunneling," JHEP 1207, 044 (2012) [arXiv: 111 1.2349 [hep-th]]. 

[5] N. Gromov and P. Vieira, "Quantum Integrability for Three-Point Functions," 
arXiv:1202.4103 [hep-th]. 

[6] I. Kostov, "Three-point function of semiclassical states at weak coupling," 
arXiv:1205.4412 [hep-th]. 

[7] R. A. Janik, P. Surowka and A. Wereszczynski, "On correlation functions of operators 
dual to classical spinning string states," JHEP 1005, 030 (2010) [arXiv:1002.4613 [hep- 
th]]. 



63 



[8] R. A. Janik and A. Wereszczynski, "Correlation functions of three heavy operators: The 
AdS contribution," JHEP 1112, 095 (2011) [arXiv: 1109.6262 [hep-th]]. 

[9] Y. Kazama and S. Komatsu, "On holographic three point functions for GKP 
strings from integrability," JHEP 1201, 110 (2012) [Erratum-ibid. 1206, 150 (2012)] 
[arXiv: 11 10.3949 [hep-th]]. 

[10] L. F. Alday, J. Maldacena, A. Sever and P. Vieira, "Y-system for Scattering Amplitudes," 
J. Phys. A A 43, 485401 (2010) [arXiv:1002.2459 [hep-th]]. 

[11] Y. Kazama and S. Komatsu, "Wave functions and correlation functions for GKP strings 
from integrability," arXiv:1205.6060 [hep-th]. 

[12] D. Gaiotto, G. W. Moore and A. Neitzke, "Wall-crossing, Hitchin Systems, and the 
WKB Approximation," arXiv:0907.3987 [hep-th]. 

[13] A. M. Polyakov, "Gauge fields and space-time," Int. J. Mod. Phys. A 17S1, 119 (2002) 
[hep-th/0110196]. 

[14] A. A. Tseytlin, "On semiclassical approximation and spinning string vertex operators 
in AdS(5) x S**5," Nucl. Phys. B 664, 247 (2003) [hep-th/0304139]. 

[15] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "A Semiclassical limit of the gauge / 
string correspondence," Nucl. Phys. B 636, 99 (2002) [hep-th/0204051]. 

[16] E. D'Hoker, D. Z. Freedman, S. D. Mathur, A. Matusis and L. Rastelli, "Extremal 
correlators in the AdS / CFT correspondence," In *Shifman, M.A. (ed.): The many 
faces of the superworld* 332-360 [hep-th/9908160]. 

[17] B. Eden, P. S. Howe, C. Schubert, E. Sokatchev and P. C. West, "Extremal correlators 
in four-dimensional SCFT," Phys. Lett. B 472, 323 (2000) [hep-th/9910150]. 

[18] M. S. Costa, J. Penedones, D. Poland and S. Rychkov, "Spinning Conformal Blocks," 
JHEP 1111, 154 (2011) [arXiv: 1109.6321 [hep-th]]. 

[19] M. S. Costa, J. Penedones, D. Poland and S. Rychkov, "Spinning Conformal Correla- 
tors," JHEP 1111, 071 (2011) [arXiv: 1107.3554 [hep-th]]. 

[20] D. Simmons-Duflin, "Projectors, Shadows, and Conformal Blocks," arXiv: 1204.3894 
[hep-th]. 

[21] F. A. Dolan and H. Osborn, "Conformal partial waves and the operator product expan- 
sion," Nucl. Phys. B 678, 491 (2004) [hep-th/0309180]. 

[22] L. F. Alday, D. Gaiotto, J. Maldacena, A. Sever and P. Vieira, "An Operator Product 
Expansion for Polygonal null Wilson Loops," JHEP 1104, 088 (2011) [arXiv:1006.2788 
[hep-th]]. 

[23] D. Gaiotto, J. Maldacena, A. Sever and P. Vieira, "Bootstrapping Null Polygon Wilson 
Loops," JHEP 1103, 092 (2011) [arXiv:1010.5009 [hep-th]]. 



64 



[24] D. Gaiotto, J. Maldacena, A. Sever and P. Vieira, "Pulling the straps of polygons," 
JHEP 1112, Oil (2011) [arXiv: 1102.0062 [hep-th]]. 

[25] L. F. Alday, B. Eden, G. P. Korchemsky, J. Maldacena and E. Sokatchev, "From corre- 
lation functions to Wilson loops," JHEP 1109, 123 (2011) [arXiv: 1007.3243 [hep-th]]. 

[26] B. Eden, G. P. Korchemsky and E. Sokatchev, "From correlation functions to scattering 
amplitudes," JHEP 1112, 002 (2011) [arXiv: 1007.3246 [hep-th]]. 

[27] B. Eden, G. P. Korchemsky and E. Sokatchev, "More on the duality correla- 
tors/amplitudes," Phys. Lett. B 709, 247 (2012) [arXiv: 1009.2488 [hep-th]]. 

[28] A. V. Belitsky, G. P. Korchemsky and E. Sokatchev, "Are scattering amplitudes dual to 
super Wilson loops?," Nucl. Phys. B 855, 333 (2012) [arXiv:1103.3008 [hep-th]]. 

[29] B. Eden, P. Heslop, G. P. Korchemsky and E. Sokatchev, "The super-correlator/super- 
amplitude duality: Part I," arXiv: 1103.3714 [hep-th]. 

[30] B. Eden, P. Heslop, G. P. Korchemsky and E. Sokatchev, "The super-correlator/super- 
amplitude duality: Part II," arXiv: 1103.4353 [hep-th]. 

[31] E. I. Buchbinder and A. A. Tseytlin, "Semiclassical correlators of three states with 
large S 5 charges in string theory in AdS 5 xS 5 " Phys. Rev. D 85, 026001 (2012) 
[arXiv:1110.5621 [hep-th]]. 



65 



